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■3 
From  high-resolution  measurement  of  melting  pressure  of  solid  -^He 

in  applied  magnetic  fields,   ranging  from  0.266   T  to  0.495   T,    the  entropy 
as  a  function  of   temperature  was  determined  by  making  use  of   the 
Clausius-Clapeyron  equation.    The  entropy  discontinuity  was  found  at   the 
transitions  between  the  paramagnetic  phase   (PP)   and  the  low-field  phase 
(LFP),    as  well   as   between  the  low-field  phase  and  the  high-field  phase 
(HFP),    indicating  that  the  nature  of   these  transitions  was  first  order. 
No  entropy  discontinuity  was  seen  at  the  transition  from  the  paramag- 
netic phase  to  the  high-field  phase,    indicating  that   it  was  not  first 
order,    but   possibly  a  second-order  transition  or   X   transition. 

The  magnetic  phase  diagram  of  solid  -%e  was  reestablished  based  on 
the  data  from  this  experiment  and  work  by  other  authors.    The 
intersection  of   the  first-order  transition  line  and  PP-HFP   transition 
line  was  determined  at  B  =  0.396   T  and  T  =  0.90  raK. 


By  fitting  the  data  of  pressure  as  a  function  of  temperature  in  the 
high-field  phase  at  B  =  0.495   T,    the  spin-wave  velocity  was  found  to  be 
6.2,    7.8,   or  8.9  cm/sec  for  one,   two  or  three  modes,   respectively.   As 
indicated  by  the  T     of  the  melting  pressure,    the  high-field  phase  was 
magnetically  ordered,   rather  than  just  a  fully  polarized  paramagnetic 
phase. 


CHAPTER   1 
INTRODUCTION 

■3 
1.1  General  Remarks  on  Solid  JHe 

Solid  -'He  is  one  of  the  quantum  crystals  (JHe,   He,  H2,  HD,  D2, 
etc.)  which  show  quantum  effects  macroscopically.  Because  of  the  small 
mass  of  He  atom  and  relatively  weak  potential  between  the  atoms 
(principally  the  Van  der  Waals  potential)  there  is  a  large  zero-point 
motion  of  solid  ^He  atoms.  The  -'He  atoms  are  not  well  localized  and 
their  wave  function  is  very  broad.  For  example,  the  amplitude  of  zero- 
point  motion  is  about  30?  of  the  nearest  neighbor  separation  at  melting 
pressure  (Roger,  Hetherington,  and  Delrieu,  1983).  When  a  first  atom 
moves  around,  second  atom  nearby  may  occupy  the  lattice  vacancy  created 
by  the  first  atom.  The  similar  process  could  happen  between  the  second 
atom  and  third  atom,  and  so  on,  causing  a  "ring"  exchange  interaction 
(McMahan,  and  Wilkins,  1975). 

On  the  other  hand,  the  nucleus  of  -'He  atom  is  a  spin  1/2  fermion. 
The  exchange  of  atom  positions  is  accompanied  by  the  exchange  of  nuclear 
spins  of  these  atoms.  The  exchange  process  in  space  among  the  JHe  atoms 
also  induces  an  effective  interaction  between  nuclear  spins.  For 
classical  solids  the  nuclear  exchange  interaction  is  always  much  smaller 

than  the  nuclear  magnetic  dipolar  interaction,  which  may  cause  a  nuclear 
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ordering  at  about  10  '  K.  But  the  exchange  interaction  in  solid  JHe 

_^ 
causes  the  nuclear  magnetic  ordering  at  10   K  (see,  for  example,  Cross 


and  Fisher,  1985).  In  fact,  solid  ^He  is  the  unique  example  which  shows 
the  nuclear  magnetic  ordering  due  to  exchange  of  whole  atom  positions. 
The  nuclear  spin  system  of  solid  ^He  provides  a  regime  with  very 
rich  magnetic  properties.  The  interaction  between  experimental  and 
theoretical  research  on  this  subject  has  much  improved  our  knowledge 
about  this  interesting  magnetic  system.  Nevertheless,  there  are  still 
ambiguities  and  disagreements  in  both  of  the  experiments  and  theories. 
Furthermore,  research  on  the  magnetism  of  solid  JHe  may  be  useful  in 
understanding  other  physical  systems,  such  as  two  dimensional  Wigner 
crystal  of  electrons,  metallic  glasses,  solid  hydrogen  and  deuterium, 
and  solid  He.  The  problem  of  magnetism  of  solid  ^He  still  atracts  a  lot 
of  attention  from  experimentalists  and  theorists. 

1.2  Heisenberg  Nearest  Neighbor  Model  and  Its  Failure 

•3 
There  are  several  basic  interaction  energies  in  solid  JHe.  The 

5 
excitation  energy  of  the  electron  from  the  closed  shell  is  of  order  10 

K.  The  characteristic  energy  of  the  phonon  interaction  in  the  solid  is 

about  10  K.  The  nuclear  dipole  interaction  has  an  order  of  1 0   K,  much 

lower  than  the  temperature  which  the  experiments  could  reach.  The  most 

interesting  interaction  in  solid  ^He,  then,  is  the  exchange  interaction. 

The  Hamiltonian  of  solid  ^He  in  an  applied  magnetic  field  could  be 

written  as 


H  =  H  ,_  +  H   +H+H.  ,  (1.1) 

ph    ex    z    d  ' 


where  H  h  is  the  part  for  phonon  interaction,  Hgx  for  exchange 
interaction,  Hz  for  Zeeman  interaction  and  Hd  for  magnetic  dipole 


interaction.  In  the  temperature  range  of  interest  the  Hamiltonian  can  be 
simplified  to  the  spin  Hamiltonian  H 

Hs  "  Hex  +  Hz  •  d.2) 

Bernardes  and  Primakoff  (1959)  first  showed  that  the  exchange 
interaction  in  solid  3He  was  many  times  larger  than  the  nuclear  dipole 
interaction.  But  the  prediction  of  the  nuclear  ordering  temperature  was 
near  0.1  K,  much  higher  than  what  the  later  experiments  indicated. 

The  obvious  theory  of  the  exchange  process  was  the  Heisenberg 
nearest  neighbor  model.  Simply  assuming  that  the  exchange  energy  J  is 
the  same  for  all  the  nearest  neighbors,  we  can  write  the  exchange 
Hamiltonian  as 

H   =  -2  J  E  I. •  I.   ,  (1 .3) 

ex  i   j   '  wo' 

i<J      J 

where  Ii  and  Ij  are  the  nuclear  spins  for  3He  atoms  labelled  by  i  and  j, 
respectively.  The  summation  is  for  all  the  nearest  pairs  in  solid  3He.  A 
negative  J  will  give  an  antiferromagnetic  structure  (Ii  and  I,  anti- 
parallel)  in  order  to  minimize  the  ground  state  energy  of  the  system.  A 
positive  J  will  give  a  ferromagnetic  structure  (Ij  and  Ii  parallel)  for 
the  same  reason. 

The  Zeeman  Hamiltonian  in  an  applied  magnetic  field  can  be  written 
as 


H  =  -  E  y  •  B 

z     .  *i 

1 


=  -  Yfi  E  S  •  B   ,  (1.4) 

i 
where  S1  and  \i     are  the  spin  and  magnetic  moment  of  the  atom  labelled  by 

i,  B  is  the  applied  magnetic  field,  and  Y  is  the  gyromagnetic  ratio. 

The  partition  function  of  solid  3He  in  an  applied  magnetic  field  is 


Z  =  Tr  exp  [  -0  (  Hex  +  Hz)  ]   ,  (1.5) 


where  8  is  1/kgT,  and  kg  stands  for  Boltzmann  constant.  The  free  energy 
can  be  obtained  by  expanding  the  partition  function  in  power  series 
of  6  (Roger,  Hetherington  and  Delrieu,  1983) 


F  -  -  r  InZ 


2       3 
?  1  [ln2  ♦  e2  f  -  e3  jL  ♦ 


\   i^f^]      d  -  66  +  c24+  '••  )  +  ..  }   .      (1.6) 


The  parameters  in  the  expansion,  e2>  e  etc.,  are  related  to  the 
exchange  energy  J.  This  expansion  is  under  the  assumption  that  the 
exchange  rate  is  small  compared  to  k  T/fi  .  The  entropy  is  given  by 


v  3T  'v 
=  NkB  (ln2  -  e2  62  /8  +  e  g3  /1 2  +  •••  )   .       (1.7) 


The  specific  heat  is 


C   =  T  (  —  ) 


Nk 

-jp  (i2  B2  -  i  B3  ♦  •••  )   •  (1.8) 


The  susceptibility  is 


P0  kB  .  92F  , 

x  =  -—7-  <  -7  ) 

3B   T 


T  -  9  +  A/T 


(1.9) 


where  the  Curie  constant  C  =  (   -5—  )      ,    yn  the  permeability  of  free 

B  Y   ft 

space,    R  the  gas  constant,   v  the  molar   volume,   and  irr—  the  Zeeman 

^kB 

splitting  in  mK  per  tesla.  The  pressure  of  solid  -*He  in  a  magnetic  field 


v  8v  ;T 
-  (ig  B/8  -  e'  B2/24  +...)+ 

+  (YhBB/2)2  (e'/2  +  a2'e/l6  +  •••)   ,  (1.10) 


where  the  primes  denote  the  derivatives  with  respect  to  molar  volume  v. 

The  first  reliable  measurement  of  the  coefficient  ep  in  zero 

magnetic  field  and  high  temperature  range  13  <  T  <  100  mK  was  done  by 

Panczyk  and  Adams  (1970).  From  the  pressure  measurement  they  found 

accurate  values  of  8e„/3v.  After  the  integration  with  respect  to  the 

molar  volume,  ~  was  determined  as  5.14  (v/24)3°.26  (mK)2  where  v  is  i 
e2 

cm^/mole.  For  the  melting  pressure  at  this  temperature  range  v  =  24.2 
cm^/mole,  that  gives  i  =  6.95  (mK)2. 


Greywall  (1977)  and  Hebral  et  al.  (1979)  determined  L  by  usi 


.ng 


specific  heat  measurements  at  high  temperature  range  20  <  T  <  50  mK.   The 
value  of  e2  from  these  experiments  was  in  good  agreement  with  the  value 
obtained  by  Panczyk  and  Adams   (1970). 

Kirk,    Osgood  and  Garber    (1969)   performed  the  first  experiment  to 
yield  the  sign  of  the  exchange  energy  by  measuring  the  magnetic 
susceptibility.    The  results  were  consistent  with  the  Curie-Weiss  law, 
giving  a  negative  value  of  8  for  Equation   (1.9),    and  thus  a  negative 
value  of  J.    Also,    the  exchange  energy  J  obtained  from  this  experiment 
was  in  good  agreement  with  the  result  of  Panczyk  and  Adams   (1970). 

In  the  period  before  1975   there  were  many  theoretical   papers  based 
on  the  Heisenberg  Nearest  Neighbor  Antiferromagnetic  model   (HNNA),   which 
gave  a  quite  satisfactory  agreement  with  the  experiments  performed  in 
the  high-temperature  range. 

The  first  important  experiment  which  contradicted  the  prediction  of 
HNNA  model  was  the   pressure  measurement  of  solid  ^He  in  an  applied 
magnetic  field  by  Kirk  and  Adams   (1971).    The  result  gave  J  =   -  0.392  mK 
at  molar   volume  23-3^   cnP/mole,    which   qualitatively  confirmed  the 
antiferrromagnetic  behavior  of  solid  ^He  with  a  negative  exchange 
energy,    but   quantitatively  only  produced  about  half  of  the  field 
dependence  predicted  by  the  theory  using  the  value  J  obtained  from  the 
early  experiment   (Figure   1). 

Another  breakdown  of  the  HNNA  model   was   the  experiment  of   the 
susceptibility  measurement.    According  to  HNNA  model   the  magnetization  M 
of  solid  ^He  would  follow  Curie-Weiss  law  as  a  first-order  approximation 


70  kG 


Figure   1.1    Pressure  measurement   in  high-field  phase 
(After  Kirk  and  Adams,    1971) 


Pressure  difference   versus   1 /T  for  molar   volume  v  =  23.34  cm3/moie 
in  different  magnetic  fields.    The  dashed  lines  are  the  theoretical 
behavior  based  on  HNNA  model.    The  solid  lines  are  the  fits  for  the 
experimental   data.    (Note:    H  =  B,    1 0  kG  =  1    T,    1    atm.    =  101.3  kPa.) 


M  =   t~TZ     •  (LID 


Prewitt  and  Goodkind   (1977)   measured  the  static  nuclear  magnetization  of 
solid  JHe  as  a  function  of  temperature  and  molar   volume  through  its 
nuclear  ordering  temperature  by  using  the  SQUID  technique.    They  reported 
that  the  magnetization  below  5  mK  increased  with  decreasing  temperature 
more  rapidly  than  the  Curie-Weiss  behavior   displayed  at  higher 
temperatures.    Another  result  found  by  these  authors  was   that   the 
magnetization  decreased  rapidly  to  40$  of   its  maximum  value  and  became 
temperature  independent   at   the  lower  temperatures.    Bernier  and  Delrieu 
(1977)   used  pulsed  NMR  technique  and  obtained  similar  results  to  those 
of  Prewitt  and  Goodkind   (1977)   above  the  transition  temperature. 

Thouless   (1965)   first  pointed  out  the  important  role  of  multiple 
exchange  mechanism  for   understanding  the  magnetic   properties  of  solid 

"3 

JHe.    He  settled  the  question  of  the  sign  for   different  exchanges,    that 
is  antiferromagnetic  for  nearest  neighbor  exchange,    ferromagnef ic  for 
three-particle  ring  exchange.    He  also  gave  the  form  for  four-spin  ring 
exchange. 

Guyer  and  his  collaborators    (1969-1975)   published  a  series  of 
papers  to   discuss  the  three-spin  and  four-spin  exchange  process    (Guyer 
and  Zane,    1969;   Zane,    1972;   McMahan  and  Guyer,    1973;   Mullin,    1975;   and 
Guyer,   Mallin  and  McMahan,    1975).    By  1974  Guyer  was  convinced  that   the 
HNNA  model  was   not  adequate  to  describe  the  magetism  of  solid  ^He. 

1.3     Magnetic  Ordering  and  the  Magnetic   Phase  Diagram 
In  the  past  decade  the  experimental  technique  for  obtaining  the 
ultra-low  temperature  below  1    mK  has  been  greatly  improved.    That  has 


provided  experimentalists  a  possibility  to  probe  directly  the  nuclear 
magnetic  ordering  of  solid  -'He. 

Halperin  et  al.    (1974)   first  observed  that  the  entropy  of  solid  ^He 
along  the  melting  curve  decreased  by  80$   in  an  interval  of  100  nK  at 
temperature  of  1.17  mK,   suggesting  a  first-order  phase  transition.   This 
result  also  demonstrated  the  failure  of  HNNA  model  which  predicted  that 
a  second-order  transition  of  nuclear  ordering  would  happen  at  about 
2  mK. 

Kummer   et  al.    (1975)    first  reported  the  determination  of  the 
entropy  of  solid  JHe  in  an  applied  magnetic  field,    and  therefore  a 
portion  of   the  magnetic  phase  diagram  of  solid  JHe.    They  found  that  at 
low  fields  below  0. 41    T  the  ordering  occured  over  a  very  narrow  tempe- 
rature interval,    and  the  ordering  temperature  was   depressed  by  the 
field.    But  above  0.41    T  the  character  of  the  ordering  suddenly  changed 
with  reduction  in  entropy  occurring  more  gradually,   and  the  ordering 
temperature  increased  with  increasing  applied  magnetic  fields.    These 
authors   interpreted  the  low-field  ordering  representing  a  magnetic  phase 
transition.   They  suggested  that  the  ordered  phase  was  probably  an  anti- 
ferromagnetic  or  spin-flop  phase.    For  the  ordering  in  higher  fields   they 
interpreted  as   paramagnetic  ordering  of  spins  by  the  applied  field, 
rather  than  a  phase  transition. 

Prewitt  and  Goodkind   (1980)    repeated  their  susceptibility 
measurements    (1977)    in  magnetic  fields  up  to  0.58   T.    Below  0.4   T  they 
obtained  the  results  more  or   less   the  same  as   that   in  the  experiment  of 
1977,    showing  that   the  ordering  was  first  order  and  that   the  ordered 
state  was  antif erromagnetic.    At  the  high  fields  above  0.4   T  they  found 
that   1/x  dropped  below  the  low-field  curve   at  a  well  defined 
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temperature,  with  the  temperature  departure  rising  when  the  applied 
magnetic  field  was  increased.  The  authors  considered  that  the  behavior 
of  1/x  reflected  the  second-order  phase  transition  between  the 
paramagnetic  phase  and  the  high-field  phase.  In  the  high-field  phase  the 
magnetization  approached  to  a  saturated  value  which  was  much  higher  than 
the  magnetization  in  the  low-field  phase. 

Nothing  about  the  magnetic  structure  of  the  high-field  phase  and 
the  low-field  phase  was  known  until  the  NMR  experiments  were  performed 
by  Adams,  Schubert,  Haas  and  Bakalyar  (1980),  and  Osheroff,  Cross  and 
Fisher  (1980)  at  almost  the  same  time  but  independently. 

Adams  et  al .  (1980)  used  a  Pomeranchuck  cell  and  worked  in  fields 
up  high  to  2.98  T.  The  %e  sample  in  the  NMR  coil  was  a  polycrystalline 
mixture  of  ordered  and  disordered  solid  at  or  near  the  transition 
temperature.  In  the  ordered  phase  below  0.41  T  they  found  resonant  peaks 
in  the  spectrum  below  and  above  the  Larmor  frequency  YB_  .  The  magnitude 
of  the  frequency  shifts  indicated  a  large  anisotropy  energy,  which  was 
incompatible  with  the  cubic  magnetic  lattice  of  most  theoretical  models 
available  at  that  time  (see,  for  example,  Hetherington  and  Willard  1975; 
Okada  and  Ishikawa  1978;  or  Roger  et  al . ,  1977).  For  the  high  fields 
above  0.41  T  Adams  et  al .  investigated  the  behavior  of  the  melting 
pressure,  magnetization  and  frequency  shift  as  functions  of  time.  From 
0.43  T  to  2.98  T  Adams  et  al.  observed  a  frequency  shift  of  about  3  kHz 
independent  of  the  fields.  This  implied  a  constant  magnetization  in 
high-field  phase  equal  to  0.55  Msat,  where  Msat  was  the  saturation 
magnetization  of  solid  %e.  The  magnetic  structure  could  not  be 
determined  by  these  data  alone.  But  from  the  result  of  this  experiment 
the  high-field  phase  was  quite  close  to  that  predicted  by  Roger,  Delrieu 
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and  Hetherington  as  a  "spontaneously  spin-flopped  state,"  which  has  two 
simple  cubic  sublattices  with  orthogonal  magnetizations.    Adams  et  al. 
(1980)   concluded  from  the  rapid  onset  of   the  frequency  shift  that  a 
second-order  phase  transition  was  occurring  at  the  high-field  phase. 

Osheroff  et  al.    (1980),   on  the  other  hand,   were  able  to  grow  single 
crystals  in  the  ordered  phase  and  study  their   properties  over  a  broad 
temperature  range  below  the  ordering  temperature.    In  principle,    NMR 
experiments  probe  the  energy  dependence  of  spatially  uniform  spin 
rotations.    The  authors  observed  large  frequency  shifts  above  the  Larmor 
frequency  at  the  low-field  phase,    with  a   zero-field  resonant  frequency 
near   zero  temperature  of  ^n/2ir  =  825  kHz.    The  large  frequency  shifts 
implied  a  large  anisotropy  of  dipole  energy  in  the  low-field  phase,   and 
a  breakdown  of  the  cubic  symmetry  in  the  magnetic  structure  of  solid 
■%e,    in  agreement  with  the  results  reported  by  Adams  et  al.    (1980). 

Furthermore,    Osheroff  et  al.    analyzed  the  spectrum  of  their 
antiferromagnetic  resonance  and  determined  the  stringent   constraints  on 
possible  sublattice  structures.    By  using  single  crystals  which  could 
support  only  three  sets  of  resonsnces,    Osheroff  et  al .    found  three 
domains   with  only  three  possible  spin  orientations.    The   direction  of 
anisotropy  was  along  one  of  the  three  principal  axes   [001],    [010], 
[100].    There  were  two  modes  observed   in  each  domain,    coupled  with  the 

existence  of   a  single  non-zero  resonance   at   zero  field.    This   implied 

2 
that   the  dipole  anisotrpy  had  the  form  E     a   (l«d)      ,   where   1  and  d  were 

anisotropy  axes   in  real   space  and  spin  space    (Figure  1.2). 

The   quasihydrodynamic  equation  of  motion  for  spin  can  be  written  as 

(Halperin  and  Saslow,    1977) 
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Figure  1.2       The   u2d2  structure 


Magnetically,    the  u2d2  structure  has   tetragonal   symmetry,   d  and 
1  are  the  spin  and  space  anisotropy  axes,   respectively. 
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S  =  YS  x  B  +  3E  /3n 


n  =  YB  -  Y   X   S   ,  (1.12) 


where  S  is  spin,  B  the  applied  field,  n  change  in  orientation  of  the 
spin  ordering  in  terms  of  small  rotations  about  three  coordinate 
axes,  ED  the  dipole  enegy  on  these  rotations,  and  ^  the  susceptibility 
tensor.  The  first  equation  is  the  torque  equation.  Notice  that  the 
dependence  of  the  dipole  energy  on  small  rotations  E  (n)  puts  an 
additional  torque  ZE^/dr]   on  the  total  spin.  The  second  equation  is  the 
kinematic  relation  between  angular  velocity  r\   and  angular  momentum.  The 
total  energy  of  spin  in  a  field  is 


1   2   ~"1 
E  =  E  (n)  +  j  Y   S  x   S  -  S-B   .  (1.13) 


The  solution  of  equation  (1.12)  under  the  condition  of  minimizing  the 
total  energy  E  was  found  by  Osheroff  et  al.  (1980)  as 


2 

\  {  WJ  +  a*  ±  [<«£-  n2Q)    *  4  u2  ^  C032Q]i/2  |  >  {un) 


where  ojl  specifies  the  Larmor  frequency,  n     the  antiferromagnetic 
resonance  frequency  at  zero  field,  9  the  angle  between  the  magnetic 
field  and  the  anisotropic  axis  in  solid  ^He  associated  with  the 
sublattice  structure.  The  resonant  spectrum  obtained  by  Osheroff  et  al. 
(1980)  was  in  good  agreement  with  the  equation  (1.14). 
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The  classic  antiferromagnetic  models  of  bcc  lattice,   such  as 
NAF,   SCAF|,   SCAF_[     ,    could  not   produce  frequency  shift  because  these 
structures  have  a  magnetic  cubic  symmetry.    The  results  of  Osheroff  et 
al.   show  that  the  original  magnetic  cubic  symmetry  for  the  low-field 
phase  was  broken.   Osheroff  et  al .    (1980)   proposed  a  possible  magnetic 
structure  for  the  low-field  phase,   u2d2  structure,   which  has  a 
tetragonal  magnetic  symmetry   (Figure  1.2).    In  the  u2d2  structure  the 
spin  orientation  alternates  to  the  opposite  direction  every  two 
successive  planes. 

Theoretical   calculation  for   u2d2  structure   (Cross,    1982;   Cross  and 
Fisher,    1985)   gave  the  value  Si  /2tt  of  1230  ty  kHz,    where  factor  ty  was   the 
renormalization  factor  due  to  zero  point  spin  fluctuations.   Compared  to 
the  experimental   value  «  /2tt  of  825   kHz,    factor  i>  was  about  0.67,    quite 
bit  smaller  than  the  value  0.85   that  derived  from  the  spin-wave 
calculation  on  simple  antiferromagnetic  states   (Anderson,    1952;    Kubo, 
1952).   More  complicated  structure  umdm   (m  is  an  integer  bigger   than  2) 
could   be  also   possible  for  the  results   of  NMR  experiments,    but   the 
longer  sequences   umdm  made  the  agreement   progressively  worse.   Whether 
the  spin  structure  is  u2d2  or   umdm  is  still  not  certain.    And  the 
quantitative  agreement   between  the  experiment  and  theory  for  u2d2 
structure  is  also  an  open  problem.    Nevertheless,    the  NMR  experiments 
have  greatly  improved  our   knowledge  about   the  magnetic  structure  of 
solid  -*He. 


*  NAF:   Normal  antiferromagnetic  phase  with  two  simple  cubic 
ferromagnetic  sublattices  with  opposite  magnetization. 

SCAF||:   Simple  cubic  antiferromagnetic  phase  with  two  simple  cubic 
antiferromagnetic  sublattices   with   parallel  magnetization. 

SCAFjj   Simple  cubic  antiferromagnetic  phase  with  two  simple  cubic 
antiferromagnetic  sublattices  with  orthogonal  magnetization. 
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Figure  1.3       Magnetic  phase  diagram  of  solid  -%e   (Experimental) 

PP:    paramagnetic   phase,    LFP:    low-field   phase,    HFP:    high-field 
phase. 

The  dashed  line  is  from  the  NMR  results  of  Adams  et  al.    The  solid 
line  is   a  fit  to  the  NMR  results  of  Osheroff  et   al.    The  dot-dashed  line 
is  the  trasition  line  determined  by  Prewitt  and  Goodkind  from  the 
measurements   of  static  susceptibility.    The  possible  magnetic  structures 
for   low-field  phase  and  high-field  phase  are  also  shown  in  the  figure. 
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Figure  1.3  is  the  magnetic  phase  diagram  of  solid  3He  based  on  the 
work  of  Adams  et  al.  (1980),  Osheroff  (1982),  and  Prewitt  and  Goodkind 
(1980).  All  the  authors  agreed  that  the  transition  between  the 
paramagnetic  phase  and  LFP,  or  between  LFP  and  HFP  were  first  order.  But 
the  order  of  the  transition  between  the  paramagnetic  phase  and  HFP  was 
not  clear.  Osheroff  et  al.  considerd  that  transition  should  be  first 
order  too.  Adams  et  al.  and  Prewitt  et  al.  found  that  was  a  second-order 
transition. 

1.4  Multiple  Exchange  Models  and  Their  Consequences 

1.4.1  The  Symmetry  Property  of  the  Multiple  Exchange  Processes 

After  the  experiments  of  Kummer  et  al.  (1975)  and  of  Halperin  et 
al.  (1974,  1978),  it  was  believed  that  the  HNNA  model  was  not 
appropriate.  There  were  several  models  proposed,  such  as  spin  glass 
model  (Beal-Monod,  1977),  vacancy  models  (Sokoloff  and  Widom,  1975a, 
1975b;  Andreev,  Marchenko  and  Meierovich,  1977;  and  Heriter  and  Lederev, 
1977,  1978),  multiple  exchange  models  (Hetherington  and  Willard,  1975) 
etc.  Some  of  these  models  are  proved  not  suitable  to  solid  ^He.  The 
multiple  exchange  phenomenological  models  were  a  success,  and  the  only 
ones  which  can  account  for  most  of  the  experimental  results  at  the 
present  time. 

The  exchange  rate  J  is  much  smaller  than  the  zero-point  motion 
frequendy  u>0  (the  attempt  frequency).  The  typical  value  of  J/u  is  -lO'1* 
from  experiment.  Therefore,  each  exchange  process  happens  rather 
separately,  and  the  total  exchange  Hamiltonian  can  be  written  as  the  sum 
of  individual  exchange  terms 
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Hov     '     £     Jn     P„  (1-15) 

(Pi 


where  Pp   is  a  permutation  operator  for  a  transposition  p,    J     is  the 
exchange  frequency  corresponding  to  P   . 

For  a  fermion  system  like  solid  ^He  the  wave  function  has   to  be 
antisymmetric  under  interchange  of  two  ^He  atoms   (both  of  spin  and 
spatial   coordinates).    Generally,   we  have 


P(R)   P(0)   =   (-DP  (1.16) 


where  P  is  the  permutation  operator  for  spatial  coordinates,  P^  is 
the  permutation  operator  for  spins,  (-1)p  is  negative  for  a  permutation 
involving  an  odd  number  of  interchanges  (even  atoms),  (-1 )p  is  positive 
for  a  permutation  involving  an  even  number  of   interchanges    (odd  atoms). 

Suppose  there  are  only  even  permutations  Pp;    then  the  total  wave 
function  i|i  must  satisfy 


Pe   *  =   +   *      •  (1.17) 


Since  the  orbital  wave  function  of   the  ground  state  ty  r  is  nodeless  from 
the  general  argument  of  quantum  mechanics    (Courant  and  Hilbert   1937), 
then  \pQ     must   be  symmetric  with  respect  to  all   allowed   permutations 


piR)  *0r-**or  •  <'-'8' 


From  Equations   (1.17),    (1.18)   we  have 
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pe°'    *0°-**0°  <'•"> 


where  xn     is  the  spin  wave  function  at  the  ground  state,   and  a  symmetric 
function  corresponding  to  a  ferromagnetic  state. 

By  a  similar  argument,   odd  permutations  will  give  an  antisymmetric 
spin  wave  function,   which  corresponds  to  an  antiferromagnetic  state. 

The  general   conclusion,   as  shown  by  Thouless   (1965),    is  that  even 
permutations  lead  to  f erromagnetism  and  odd  permutations  lead  to 
antiferromagnetism.    This  symmetry  property  is  an  important   criterion  for 
determining  what  kind  of  multi-exchange  interaction  takes  place  in 
different  magnetic  phases  of  solid  JHe. 

1.4.2  The  different  exchange   processes  and  their  consequence  for   bcc 

solid  -*He 

The  Lennard-Jones   interaction  between  -'He  atoms  shows  a  very  deep 

repulsive   potential  and  very  weak  attractive   potential.    The  JHe  atom, 

therefore,    could  be  treated  as  a  hard  core  sphere  with  diameter 

au  =  2.14  A.    The  hard   core  may  cause  two  spin  exchange  to  have  a 
ri 

tremendous   large  effective  potential   compared  to  ring  exchange  of 
nearest  or  even  next  nearest  neighbors.    That  is  the  basic  physical 
argument  of  the  multiple  exchange  models.   We  consider  only  the  most 
compact  3-particle  and  4-particle  cyclic  exchange  in  the  following 
analysis. 


Three-particle  cyclic  exchange.    Any  permutation  operator  may  be 
expressed  as   the   product  of  pairwise   permutation,   which  in  turn  can  be 
written  in  terms  of  Pauli  matrices  as 
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Pij    =   2    (1    +   V    °j)       •  ^1-20) 


The  3-particle  exchange  process  is  shown  in  Figure  1.4a.   The 
operator  for  3-particle  cyclic  exchange   (i  +  j   ->  k)    is,    then 


P.  .,    =   P.  .    P., 
ijk         ij     jk 

=  4    0    +  o. •    a.)    (1    +  o.«    a J      .  (1.21) 


i       J  j        k' 


By  using  the  relation 


(V    °j)    (V    V   =  °j*    °k+   i  V    (oj   x  V  ^-22) 


we  have 


pijk  -U' +  v  °j +  v  °k  +  v  °i + 


There  is  also  a  conjugate  exchange    (i  +■  j   <-  k)    ,   which  has  the  fc 


Pijk"1   -  H1    +  V   W   °k  +  V   °i  - 

-    i   o.-    (a.    x   ok)]      .  (1.24) 


The  sum  of  Equations    (1.23)    and   (1.24)    gives   the  3-particle  cyclic 
exchange  operator   as 


P.  .,     +   P.  .,     1 
ijk  ijk 

=  2    (1    +   V    °j    +   °j*    °k  +   V    °i)      •  C1.25) 
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The  Hamiltonian  for  3-particle  exchange  is 


T 
Hex   -   Jt        *        ^i^'h,      )  d.26) 


eX         fc  <ijk>        ^         «* 


which  could  be  reduced  to  the  sum  of  pairwise  spin  interaction  according 
to  Equation  (1.25).   So,   the  3-particle  exchange  alone  leads  to  the  same 
universality  class  as   the  pairwise  interaction   (i.e.   Heisenberg 
interaction).   We  would  expect  that  3-particle  exchange  process  may  cause 
only  a  second-order  phase  transition  the  same  as   the  Heisenberg  nearest 
neighbor   interaction.    There  is  no  doubt  about   the  existence  of  the 
first-order   phase  transition  between  the  paramagnetic  phase  and  the  low- 
field  phase,    as  all   the  experimental  results  have  shown.    Also,    the  3- 
particle  exchange  process  favors  f erromagnetism  on  contrast  to  the 
antiferromagnetic   property  of  the  low-field  phase. 

Four-particle  cyclic  exchange.    The  4-particle  cyclic  permutation 
operator  Pijkl   can  be  written  as 

Pijkl  -   Pijk   *    Pn      •  (1-27) 

By  using  Equations   (1.20)   -   (1.25),    the  4-particle  exchange  operator  can 
be  expressed  as 


P  +  P       _1 

ljkl         ijkl 

■  jT   (1    +        I       o-    a.   +   G.  .,  ,  )  (1.28) 

4  .  i        ]  l ikl 

pairs  J  J 


21 


where   E   is  taken  over  the  six  distinct  couples  among  the  four 
pairs 

particles  (i,  j,  k,  1),  and 


Gijki  -  (v  V  (v  V  + 

+  (cy    ax)    (a  •    ok)   -   (o^    ok)    (a  •    o1)      .  (1.29) 


The  operator  G^^  can  not  be  reduced  further  to  Heisenberg  pairwise 
operator.    It  is  this  operator  which  leads  to  a  first-order  phase 
transition. 

Roger,   Hetherington  and  Delrieu  (1983)   studied  the  4-particle 
exchange  processes  in  great  detail  by  using  a  computer  to  minimize  the 
free  energy  of  the  system.    They  calculated  two  kinds  of  most  compact   4- 
particle  cyclic  exchange   process  involving  four  first  neighbors.    One  is 
folded  4-particle  exchange,   the  other  is  planar  4-particle  exchange 
(Figure   1.4b).    It  turns  out  from  their   calculation  that  for  folded  4- 
particle  exchange  the  magnetic  structure  should  be  the  simple  cubic 
antiferromagnetic   (SCAFJ_,    Figure   1.5a).    Because  of  the  cubic  symmetry  of 
SCAF_|_  structure,    the  anisotropic  component  of  the  dipole-dipole 
interaction  between  the  nuclear  magnetic  moments   vanishes.    This 
contradicts   the  experimental  result  of  a  very  high  antiferromagnetic 
resonance  frequency  which  requires  a  large  dipolar  anisotropy. 
Therefore,    the  SCAF_|_  structure,    induced  by  folded  4-particle  exchange 
alone,    is  not  adequate  for  describing  the  magnetic  structure  of  the  low- 
field  phase.    For  the  planar  4-particle  exchange  Roger  et   al.    found  a 
stable  structure  with  minimum  free  energy,    simple  square  antiferromag- 
netic phase   (SSQAF,    Figure   1.5b).    Although  the  SSQAF  structure  has  a 
large  dipolar  anisotropy,   which  produces   two  degenerate  longitudinal 
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(a) 


b) 


Figure   1.4       Compact   exchanges   in  bcc  lattice 

(a)  The  most   compact  3-particle  cyclic  exchange  in  a  bcc  lattice. 

(b)  The  most  compact   4-particle  cyclic  exchange   in  a  bcc  lattice. 
(1-2-3-4)  represents  a  folded  exchange.    (1-2 '-3* -4')   represents  a  planar 
exchange. 
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(b)    SSQAF 

Figure  1.5   SCAF|_  and  SSQAF  structures 

(a)  Magnetic  structure  of  simple  cubic  antiferromagnetic 
phase  (SCAFj_),  which  consists  of  two  interpenetrating  simple  cubic 
antiferromagnetic  structures  with  perpendicular  magnetizations. 

(b)  Magnetic  structure  of  simple  square  antif errromagnetic  phase 
(SSQAF).  The  lines  [100]  are  ferromagnetic.  Perpendicular  to  [100]  there 
is  a  planar  structure  with  two  interpenetrating  simple  square 
antiferromagnetic  lattices. 
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antiferromagnetic  resonance  modes  and  one  transverse  mode,   the  model 
still  cannot  represent  the  observed  phase.    According  to  the  theoretical 
calculation  of  Roger  et  al.    (1983),   the  transverse  mode  has  resonance 
frequency  1.4  MHz,   much  larger  than  the  observed  value  of  825  kHz  in  the 
experiment  of  Osheroff  et  al.    (1980).   The  longitudinal  modes  will 
separate  into  two  modes  in  a  magnetic  field,   giving  a  more  complicated 
spectrum  than  that  observed  by  Osheroff  et  al.    (1980). 

Two-parameter  model.   Since  none  of  the  single  parameter  models  can 
explain  the  experimental  data  satisfactorily,    it  is  natural   to  think  of 
multi- parameter  models.    The  philosophy  is  to  find  the  Hamiltonian  with 
the  fewest  parameters  which  could  account  for  the  experimental  data. 

The  most  successful  model   in  the  two-parameter  model   proposed  by 
Roger  et  al.    (1980,    1983)    includes  3-particle  exchange  Jfc  and  planar  4- 
particle  exchange  K   .    The  Hamiltonian  of  this  model  can  be  written  as 

T  -1 

H        -J.         I        (P.  .,     +   P.  .,       )    - 
ex         t  <ijk>        xjk         ijk 

P  -! 

"  KP         l       (PiJkl   +  Piikl      J  (1'30) 

p  <ijkl>       1JK1         1Jkl 

where  the  sum  is  taken  over  all   distinct   3  or   4-particle  exchanges. 

According  to  symmetry  property  of  the   ground  state,    Jf   and  Kn  both  are 

z  p 

negative.    The  structure  was  solved  by  computer  simulation.    The  ground 

state  is  an  up-up-down-down  phase    (u2d2)    under  the   condition    |jJ    > 

0. 25 | Kp | ,    in  agreement  with  the  model   propsed  by  Osheroff  et  al.    (1980). 

Roger  et  al.    (1983)   evaluated  the  values   of  Jh   and  K„  by  using  the 

t.  p 

coefficient  e2  and  spin-wave  velocity  v,   which  can  be   calculated  exactly 
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with  4-spin  exchange.   Compared  with  the  experimental  values  of  e     and  v, 

the  Jt  was  found  to  be  -0.13  mK,   and  K     -0.385  mK. 

The  other  properties  could  be  calculated  from  the  set  of  Jf   and  K   . 

c     p 

Table  1.1  lists  a  comparison  between  the  results  calculated  by  Roger  et 
al.  (1983)  and  experimental  values.  They  are  in  quite  good  agreement. 

Table  1 
Comparison  of  u2d2  model  with  experimental  results 


Experimental  results     Calculated  results 

Jt=-0.13mK,  Kp=-0.385mK 


Susceptibility  c/5.9   <   x   <  c/5.2  x   =  c/5.08 

(Morii   et  al . ,    1978) 

Zero-field  resonance  ft     =  825   kHz  800   <  fl     <  860  kHz 

frequency  (Osheroff  and  Yu,    1980)  ° 

Mean  spin-wave  c  =  8.4  ±  0.4   cm/sec  8  cm/sec 

velocity  (Osheroff  and  Yu,    1980) 


Figure   1.6  shows   the  phase  diagram  based  on  two-parameter  model 
(Jf   Kp)'    There  is  a  first-order  transition  between  the  low-field  phase 
and  the  paramagnetic  phase  at  low  magnetic  fields,   and  a  first-order 
transition   between  u2d2  phase  and  the  high-field  phase  at  about    1.6   T 
below   1    mK.    Also  the  phase  diagram  shows  a  second-order  transition 
between  the  high-field  phase  and  the  paramagnetic  phase  at  high  fields. 
These  features  are  qualitatively  consistent  with  the  experimental  phase 
diagram  shown  in  Figure   1.2. 

There  is  another  first-order  line  ending  at  an  intersection  point 
and  extending  upward.    This  is   consistent   with  the  results  reported  by 
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Figure  1.6       Phase  diagram  of   two-parameter  model    (theoretical) 


The  phase  diagram  is   based  on  two-parameter  model   of  Roger   et  al 
(1983)   with  J.   *    -0.13  mK  and  K_        -O.^fis  m*.    The  solid  lines   are  the 


.t  -   -0.13  mK  and  K     =   -0.385  mK. ,  ,, 

first-order  trasition  lines.    The  dashed  line   predicts  a  second-order 


trasition.   The  dot-dashed  line  in  the   insert  is  also  a  second-order 
trasition  line   predicted  by  Heisenberg  nearest  neighbor 
ant i ferromagnetic  model.    The  proposed  magnetic  structures   for  the  low- 
field  phase  and  high-field  phase  are  also  shown  in  the  figure. 


27 


Osheroff  (1982),  although  it  is  not  clear  experimentally  how  and  whether 
this  line  ends.  The  theoretical  calculation  also  predicts  that  the 
second-order  line  bends  over  at  very  high  magnetic  fields  of  about  16  T. 
No  experiment  has  reached  such  high  fields  yet. 

There  are  some  other  models  calculated  by  Reger  et  al.  (1983)  such 
as  a  three-parameter  model  including  3-particle  exchange  (J*.),  both 
planar  and  folded  4-particle  exchanges  (Kp,  Kf);  or  a  three-parameter 
model  including  2-particle  (JNN) ,  3-particle  (Jfc)  and  planar  4-particle 
(Kp)  exchanges.  These  models  do  not  fit  the  experimental  results  as  well 
as  two-parameter  model  (Jt,K  ). 

Physically,  an  exchange  process  depends  on  the  free  space  available 
without  causing  a  drastic  potential  increase  between  particles.  A 
puzzling  question  is  why  different  exchange  processes  could  happen  with 
comparable  rates  at  the  same  molar  volume  of  solid  %e.  Unless  there  is 
an  yet  unknown  mechanism  which  can  yield  many  exchange  processes  at 
comparable  rates,  it  is  unlikely  that  several  different  exchange 
processes  will  happen  in  the  same  system  with  comparable  rates.  A 
candidate  theory  to  explain  the  different  multiple  exchanges  having 
comparable  rates  is  "zero  point  vacancies"  which  is  a  finite 
concentration  of  vacancies  existing  in  the  ground  state.  These  vacancies 
may  play  an  essential  role  in  the  exchange  mechanism  (see  discussion  of 
Cross  and  Fisher,  1985,  for  detail). 

1.5  Purpose  of  This  Work 
Most  3-dimensional  localized  magnetic  systems  show  only  small 
quantum  effects  that  can  be  accounted  for  by  the  spin-wave  perturbation 
theory  (Anderson,  1952;  Kubo,  1952).  Only  solid  %e  provides  a  highly 
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quantum  3 -dimensional  magnetic  system,   which  we  expect  to  understand 
from  first  principle  of  quantum  mechanics.    In  the  past  ten  years 
experimental  and  theoretical   physicists  made  great  efforts  to  build  up  a 
picture  of  the  magnetism  of  solid  -'He  both  in  experimental  and 
theoretical  views  and  to  explain  this  system  in  a  consistent  way.   Up  to 
today  there  are  still   questions,   uncertainties  and  disagreements  in  this 
field.    It  will  be  always  controversial  until   the  physics  about  it  is 
completely  understood  and  no  longer   attractive  to   physicists. 

Although  the  multiple  exchange  models  seem  to  account   best  for  most 
of  the  experimental  data,   there  are  still   quantitative  discrepancies 
between  the  experimental  phase  diagram  and  the  theoretical   phase 
diagram.   For  theorists,   more  accurate  calculation  beyond  the  mean  field 
approximation  is  needed  to  verify  the  correctness  of  the  multiple 
exchange  models.    In  the  experimental   phase  diagram  the  region  around  the 
intersection  of  phase  lines   is  quite  ambiguous.   The  following  are  the 
main  discrepancies  among  the  theories  and  experiments. 

1.  The  multiple  exchange  models   (2-parameter ,    3-parameter)   predict 
a  first-order  phase  line,   extending  upward  and  ending  at  an  intersection 
point.    The  theory  also  predicts  a  second-order   phase  line  adjacent  to 
the  previous  one   (first-order  line)   but  bending  over  at  high  magnetic 
fields    (see  Figure  1.6).    None  of   the  experimentl   performed  before  found 
both  of  these  two  lines.    Either   a  first-order   phase  line  or   a  second- 
order  phase  line  was  reported. 

2.  Experimentally,   the  transition  between  the  paramagnetic  phase 
and  the  high-field  phase  has   been  reported  as  first  order    (Osheroff, 
1982;   Cross,    1982),    second  order    (Prewitt  and  Goodkind,    1980;   Kummer  et 
al.,    1975;   Adams  et  al.,    1980),    or   A   transition   (Uhlig  et  al.,    1984).    It 
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raises  a  question:    Is  the  high-field  phase  a  totally  polarized 
paramagnetic  phase  or  a  new  phase?  Roger,   Hetherington  and  Delrieu 
suggested  that  the  high-field  phase  was  a  canted  normal  antiferro- 
magnetic  phase   (CNAF).   This  state  has  a  broken  symmetry  transverse  to 
the  magnetic  field  and  is  separated  from  the  paramagnetic  phase  by  a 
second-order  phase  line. 

The  purpose  of   this  work  is  to  probe  the  region  around  the 
intersection  point  of  different   phase  boundaries  in  the  phase  diagram, 
to  confirm  the  order  of  the  transition  between  the  paramagnetic  phase 
and  the  high-field  phase,   and  to  determine  whether  the  high-field  phase 
is  an  ordered  phase  or  just  paramagnetic  ordering  by  the  applied  field. 
We  hope  that  the  information  from  this  work  would  remove  the  ambiguity. 

A  sensitive  way  of  finding  the  entropy  of  solid  3He  and  then  the 
order  of  the  transition  is   by  using  the  melting  pressure  and  the 
Clausius-Clapeyron  equation 


IP  VS1  ,  , 

dT         V   -  V  '  d  -31  ) 

s        1 


where  Ss>    S1,    Vs  and  V1  are  the  solid  and  liquid  entropies  and  molar 
volumes,   respectively.   The  method  used  in  this  experiment  is  to  make 
high-resolution  measurement  on  the  melting  pressure  at   different  applied 
magnetic  fields  P(T,B)    in  the  region  0.2   <  B  <   0.5   T,    spanning  the  low- 
and  high-field  phases,   and  in  a  temperature  range  going  well   below  the 
transition  temperature.    A  specially  designed  sample   cell  serves  for 
pressure  measurement.    A  pulsed  NMR  thermoneter  measures   the  temperature 
of  the  cell. 


CHAPTER  2 
APPARATUS 

2.1  General  Description 
The  magnetic  ordering  transition  of  solid  3He  at  melting  pressure 
happens  at  about  1  mK.  A  dilution  refrigerator,  combined  with  nuclear 
demagnetization  technique,  is  able  to  cool  3He  at  melting  pressure  down 
to  below  1  mK  without  difficulty. 

In  order  to  perform  an  experiment  at  such  a  low  temperature  several 
general  experimental  principles  have  to  be  followed  (see,  for  example, 
Lounasmaa,  1974).  In  this  experiment  the  cryostat  and  most  of  the 
electronic  instruments  rest  in  a  40.5  m3  (5  m  long  x  3  m  wide  x  2.7  m 
high)  copper  screen  room  which  is  grounded  to  the  earth  at  one  single 
point.  Low  pass  filters  for  0-60  Hz,  1 20  V  AC  power  line  and  high 
frequency  filters  for  eliminating  the  burst  at  3150  Hz  used  for  campus 
clock  synchronization  are  installed  to  give  a  "clean"  power  supply  for 
all  the  instruments  inside  the  screen  room.  All  vacuum  pumps  and  the  gas 
handling  board  of  dilution  refrigerator  are  outside  the  screen  room  as 
well  as  the  computer  and  the  data  acquisition  unit.  Great  care  has  been 
taken  to  isolate  the  cryostat  from  vibration  sources.  For  example,  the 
dewar  is  suspended  from  a  triangular  aluminium  plate  which  is  supported 
at  its  three  corners  by  pneumatic  isolation  mounts,  model  XL-A  made  by 
NRC  (Newport  Research  Corporation).  Bellows  are  put  in  pumping  lines 
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between  the  cryostat  and  pumps.    These  methods  are  proved  to  be  very 

_Q 

effective  for  reducing  the  heat  leak  down  to  the  order  of  1 0   J  W. 

Figure  2.1   shows  the  gross  dimensions  of  the  dewar  and  cryostat. 
The  dilution  refrigerator,   the  PrNic  bundle  for  nuclear  demagnetization 
and  experimental  volume  are  inside  the  vacuum  can,   as  shown  in  Figure 
2.2.    The  average   consumption  of  liquid  helium  is  about  22  liters   per 
day,    including  transfer  losses. 

A  home-made  superconducting  magnet  for  the  nuclear  demagnetization 
is  made  of  NbTi  matrix  wire.    The  maximum  field  at  the  center  of  this 
magnet  can  go  up  to  5  T.    A  commercial  NbTi  superconducting  magnet,   made 
by  NALORAC   provides  a  magnetic  field  for  NMR  thermometer  and  -*He  sample 
with  high  homogeneity  of  0.3  ppm  at  the  experimental  region. 

2.2     Dilution  Refrigerator 
A  dilution  refrigerator  has   become  standard  equipment  for 
experimental   research  below   1    K  since  the   principle  of  the  dilution 
refrigeration  was  suggested  by  London   (1951). 

A  remarkbale  property  of  JHe  and     He  mixture  was   discovered  by 
Edwards  et  al.    (1965).    It  is   of  most  importance  that  the  equilibrium 
concentration  of  JHe  in  the  dilute  phase  of  JHe  and     He  mixture  happens 
to   be  finite.    Even  at  absolute   zero  temperature  this   concentration  of 
-'He  is  still  6.4?.    Below  0.5  K  liquid     He  is  effectively  in  its   quantum 
mechanical  ground  state  due  to  its   zero  nuclear  spin  and  superfluid 
properties.    There  are  practically  no  phonon  or  roton  excited  in  liquid 
^He.    Therefore  liquid     He  has  no  entropy  or  heat   capacity  for  the 
cooling  process.    In  contrast  to     He,    -"He  atom  with  nuclear  spin  1=1/2 
obeys  Fermi-Dirac  statistics.    Landau's  theory  of  Fermi  liquid  applies  to 


32 


Jin 


Vapor-cooled. 
Current  Lead 


Baffles 


7 


ss 


Vaccum  Can  Pumping  Line 


Superconducting  Magnet 
for  Nuclear  Demag. 


30  cm 
I 1 


ss 


1 


y 


\\\\\\\\\\\\y\\\\\\\*A 


Support    Table 


IK  Pot  Pumping  Line 


lHc  Pumping  Line 


Vaccum  Can 


Superconducting  Magnet 
for  NMR  and     3He 
Sample 
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liquid  JHe  in  the  concentrated  phase   (Wilkes,    1967).   The  entropy  and 
heat  capacity  of  liquid,    as  a  normal  Fermi   liquid,    is   proportional  to 
temperature  T  near   absolute  zero  degree,   providing  a  significant  cooling 
power  in  the  temperature  range  which  is  of   interest   in  dilution 
refrigeration. 

The  main  part  of  the  dilution  refrigerator  used  in  this  work  is  a 
commercial  one,   model  DRP-43  made  by  former  SHE  company.   The  3He-4He 
solution  used  in  the  refrigerator  contains   1.35  moles  of  ^He  and  3. 64 
moles  of     He,   equivalent  to  50   cm3  of  ^He  and  100   cm3  of  ^He.    Basically, 
all   the     He  is  contained  in  the  dilution  refrigerator  unit.    The  ^He  in 
the   dilution  refregirator  is  about   1.1    moles   at  a  condensing  pressure  of 
100  torr.    The  rest  of   the   ^He  is  circulating  by  the   pumping  system. 

Figure  2.2  shows   the  dilution  refrigerator   unit.    Incoming  ^He  gas 
condenses   into  liquid  at  the   1    K  pot,   which  keeps   at  about   1.3  K  by 
pumping     He  in  the  pot.   The  liquid  DEe  then  comes  into  the  still,   which 
cools  the  liquid  3He  further   down  to  0.5-0.7  K.    A  heater   inside  the 
still   is  designed  for  adjusting  the  3He  flow  and  cooling  power.    After 
passing  through  a  continuous  heat  exchanger   and  three  discrete  heat 
exchangers,    the  liquid  reaches  mixing  chamber,    the   coldest   part   in  the 
dilution  refrigerator   unit.    The  cooling  process   in  the  mixing  chamber 
happens  at   the   phase  boundary  due  to  ^He  atom  transferring  from  ^He 
concentrated  phase  into  the  dilute  phase.    The  cooling  power  Q  is 

proportional   to  the  JHe  circulation  rate  n_   and  the  square  of   the  mixing 

•         •         2 
chamber  temperature   (Q  a  n     T   ).   Table  2.1    lists  some  operating 

parameters  of  the  dilution  refrigerator  used  in  this  experiment.    The 

typical  JHe   circulation  rate  is  about   0.5  mmole/sec.    It  takes  24  to  36 


Figure  2.2   Experimental  arrangement 

A  1K  pot 

B  still 

C  continuous  heat  exchanger 

D  discrete  heat  exchangers   (three  stages) 

E  mixing  chamber 

F  Tin  heat  switch  between  the  dilution  refrigerator  and 
nuclear  demagnetization  stage 

G  vespel  supporting  rods 

H  thermal  link  (Cu)  to  the  nuclear  demagnetization  stage 

I  superconducting  magnet  for  the  nuclear  demagnetization 

J  PrNic  nuclear  demagnetization  stage 

K  superconducting  magnet  for  NMR  thermometer  and  JHe  sample 

L  thermal  link  (pure  Ag)  to  %e  sample  cell  and  strain  gauge 

M  Pt  pulsed  NMR  thermometer 

N  -^He  sample  cell 

0  vacuum  can 

P  strain  gauge 

Q  thermal  shield  linked  on  mixing  chamber 

R  thermal  shield  linked  on  still 
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hours   to  precool   the  nuclear   demagnetization  stage    (0.6  moles  of  PrNir) 
down  to  about  8  mK  in  a  magnetic  field  of  4   T. 

Table  2 
Operating  parameters  of  DRP-43 


Heat  load  of  the  mixing  chamber    (mW)  0.15  0.4 

Still   power    (mW)  20  20 

-%e  circulation  rate   (mmole/sec)  0.75  0.83 

Mixing  chamber  temperature   (mK)  20  95 


2.3     Nuclear  Demagnetization  Stage 

Magnetic  ordering  of  solid  %e  at  melting  pressure  occurs  at  about 
1    mK.    In  order  to   perform  the  experiment   at  such  a  low  temperature 
cooling  methods  other  than  the  dilution  refrigerator   are  necessary. 
Pomeranchuk  cooling  and  nuclear  demagnetization  are  the   commonly  used 
techniques   to  reach  the  temperature  range  of  millikelvin  or  sub- 
millikelvin. 

The  principle  of  nuclear  demagnetization  was  proposed  by  Gorter 
(1934).    Because  the  nuclear  magnetic  moments  are  about  2000   times 
smaller  than  the  electronic  magnetic  moments,    2000  times  larger   values 
of  B^/Ti   are  required  in  order  to  obtain  the  same  entropy  reduction  in 
the  nuclear  system  than  in  the  electronic  system,   where  B,    stands  for 
the  initial  magnetic  field  and  T,   stands  for  the  precooling  temperature. 
Nuclear  demagnetization  became  practical   only  because  of   the  succeccful 
development  of  the  dilution  refrigerator  and  superconducting  magnet. 

For  a  non-interacting  nuclear  dipole  system  in  a  magnetic  field  the 
partition  function  is 
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rm=+I  nNn 

Z  =    [    I_     exp   (yngnmB/kBT)]  ,  (2.1) 


where  I   Is  the  nuclear  spin,    un  the  Bohr  nuclear  magneton,   gn  the 
nuclear  Lande  g-f actor,  m  the  magnetic  quantum  number,   B  the  external 
magnetic  field,   n  the  number  of  moles  of  the  nuclear  specimen,    NQ 
Avogadro's  constant,   kB  Boltzmann  constant. 

The  entropy  of  the  nuclear   dipole  system,    therefore,    is 

S  -  nR  {  -f  coth^f  -  I^US  coth^glix  + 

♦  in   [   sinh^illx   /sinh_x   ]}      ^  (    2>2) 

where  x  =  U^mB/kgT  and  R  is  the   gas   constant. 

Under  the  condition  of   the  adiabatic  demagnetization  B/T  must  be  < 
constant,  which  gives 


Tf  =  bT  Ti      •  (2.3) 

1 


where  the  subscript  i  refers  to  the  initial  state  and  f  refers  to  the 
final  state. 

Considering  the  internal  dipolar  interaction,  Equation  (2.3)  should 
be  modified  to 

B  2  ♦  b2 

where  b  represents   the  effective  dipolar  field   (see  Lounasmaa,    1974). 
The  heat   capacity  and  the  magnetic  susceptibility  can  be  written  as 
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2  2 

-  +  t>  : 

B   =  3 '  (2-5) 
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,,     _      A    (B"   +  b    ) 


X  =  y         (Curie's  law)      ,  (2.6) 


where  uQ  =  free  space  permeability, 

A   =   nN0I(I+1)    M0Pn2gn2/3kB  (2.7) 

The  method  of  hyperfine  enhanced  nuclear  demagnetization  is  also 
widely  used  for  reaching  millikelvin  temperature  (Altshuler,    1966; 
Andres  and  Bucher,    1968).    A  large  hyperfine  field  can  be   induced  by 
moderate  external  fields  in  singlet  ground  state  ion  with  high  Van  Vleck 
susceptibility.    The  enhancement  factor   is  defined  by  a  =  1  +K  where  K   is 
the  Knight  shift.    The  range  of  the  factor  a   is  about  from  10  to  200, 
giving  a  large  effective  field.    All  the  formulas  of  Equation   (2.2)    to 
(2.7)   are  still   effective  with  the  replacement  of  B  by  oB. 

The  nuclear   demagnetization  stage   in  this  experiment   uses  0.6  moles 
of  PrNi5   as  refrigerant;   PrNi5  has  Pr+5   ion  which  can  produce  a 
hyperfine  enhanced  magnetic  field  under  the  interaction  of  an  external 
magnetic  field.    The  factor  a  of  PrNi5  is   16.4.    The  nuclear  ordering 
temperature  of  PrNig   is  near  0.4  tnK,   which   is  about  the  temperature 
limitation  by  nuclear  demagnetization  of  PrNi5    (Kubota  et   al.,    1980).    An 
important  difference  between  the  hyperfine  nuclear  enhahced 
demagnetization  and  the  brute  force  nuclear   demagnetization  on  copper  is 
the  much  shorter  nuclear  spin-lattice  relaxation  times   t     that  are 
encounted  in  hyperfine  enhanced  demagnetization  meterials;    x      is  usually 
so  short    (of  order   10  us  at   1    K)   that   it  can  not   easily  be  observed 
experimentally.    So  the  electrons  are  always   in  local   thermal   equilibrium 
with  the  rare  earth  nuclei. 
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It  is  worth  mentioning  that  some  meterials  other  than  copper,   PrNit- 
etc.   are  very  promising  for  being  nuclear  refrigerants.    For  example, 
indium  is  a  very  suitable  material   for  the  first  demagnetization  stage 
in  a  two  stage  demagnetization  cryostat  because  of  its  large  Curie 
constant,    small   Korringa  constant,   high  filling  factor  for  making  the 
nuclear  demagnetization  stage  and  its  mechanical   property  for  easy 
handling  (Tang  et  al.,    1985). 

2.k     Thermal   Isolation 
In  order  to  stay  at  millikelvin  temperatures   for   a  long  time  to  do 
the  experiment,    it  is   particularly  important  to   prevent  heat  flow  into 
the  experimental   region  and  to  reduce  the  heat   leak  down  to  a  level  of 
nanowatts. 

The  section  2.1    has   described  several   effective  ways  to  cut   down 
the  heat  leak  from  the  mechanical   vibration  and  electromagnetic 
radiation.    The  other   principal  source  of   the  heat   leak  is  from  the 
thermal  conduction  through  supports  or  residual   gas  in  the  vacuum  can. 
It  is  a  common  practice  to  anchor  heat  shields,   supports,   tubes, 
electric  leads  to  intermediate  cooling  stages.    In  our   cryostat  two  heat 
shields  are  anchored  to  the  still   and  the  mixing  chamber,   respectively. 
These  heat  shields  are  made  of  0.08  cm  thick  oxygen  free  high  conductive 
copper    (OFHC  copper).    The   bottoms  of  the  shields  are  aligned  with  each 
other   by  a  series   of  spacers  made  of   brass  ring  and  thin   vespel   pieces 
as   isolated  parts.    These  two  heat  shields  effectively  absorb  the 
blackbody  radiation  from  the  H  K  helium  bath. 

Tubes  reaching  the  low  temperature  region  for   ^He  sample  and 
coaxial   cables  are  made  of   copper-nickel  which  has   poor  thermal 
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conductivity  and  quite  good  electrical  conductivity  at  low  temperatures. 
These  tubes  are  anchored  at  all  the  different  temperature  stages  from 
1    K  pot,   still,   first  step  heat  exchanger  to  mixing  chamber.   Advance 
wire   (copper-nickel  alloy)   is  used  for  the  electrical  leads.    For  coaxial 
cables  Apiezon  N  grease  is  introduced  into  the  tubes   to  improve  the 
thermal  contact  between  leads  and  tubes  in  order  to  reduce  the  heat  leak 
along  the  leads. 

The  residual  exchange  gas   in  the  vacuum  can  is  possibly  a  serious 
source  of  the  heat  leak.    It  has  been  reported  that  the  ortho-para 
conversion  of  hydrogen  molecules  causes  a  time  dependent  heat  leak  up  to 
100  nW   (see,   for  example,    Berg,    1983).    This  heat  leak  will  decrease  to 
an  acceptable  level  after  a  couple  of  weeks  from  the  start  of  cooldown. 
In  order  to  avoid  this  ortho-para  conversion  of  hydroden  we   use 
1000  ytorr     He  gas  as   the  exchange  gas,   which  is  pumped  out  for  about  6 
hours  at  temperature  of   10  K.    After  a  couple  of  nuclear  demagnetization 
the  heat  leak  stabilizes  at   1    nW. 

2.5     Gas  Handling  System 
Fig-are  2.3  shows   the  gas  handling  system  for  the  dilution 
refrigerator.    The   general  feature  of  the  system  could   be  referred  to 
Lounasmaa's  book   (1974).   We  only  mention  several   special   poinsts  for  our 
cryostat. 

The  filter  F1    in  the  path  of  9B3   boost  pump  is  designed  for 
preventing  oil   vapor   of  9B3  and  other  materials   of  the   cracked  oil   from 
backstreaming  into  the  3He  circulation  sy3tem>    The  flUep   ig  a 
cylindrical   container   (9.5  cm  in  diameter,    15  cm  high)   full  of  molecular 
sieves    (Linde  Adsorbents,    Type   13X,    1/16  pellets).    In  order  to  prevent 


42 


—  o 

<> 

ff 

a  v 

ZD 

E 

3 

a. 

o 

c 

" 

-si 

43 


oil  cracking  the  heater  power  of  9B3   is  reduced  to  60$  of  the  recom- 
mended value.   The  liquid  N2  cooled  traps  F3  and  F4,   located  in  the  3He 
condensed  line  after  the  mechanical  pump  ED660,   have  a  similar  function 
as  filter  F1 .   Only  one  of  the  F3  and  F4   is  in  use.    The  other  is  ready 
for  replacement  when  one  is  going  to  be  blocked.   Uhlig  et  al.    (1983)   has 
described  a  special  design  of  the  traps  F3  and  F4.   With  the  help  of 
these  measures   the  cryostat  has   been  operated  properly  for   10  monthes 
without  blocking. 

Figure  2.4   is  the  diagram  of  3He  gas  handling  system.    By  using  the 
dipstick  along  with  liquid  4He  storage  dewar,    it  is  convenient  to 
introduce  3He  from  the  tank  at  room  temperature  into  the  sample  cell  and 
strain  gauge  in  the  cryostat.    The  dipstick  is  made  of  a  stainless  steel 
tube   (1.27  cm  in  diameter,   20  cm  long)   filled  with  charcoal.    It  can 
easily  raise  the  pressure  in  the  cells  up  to  10  MPa. 

2.6     JHe  Sample  Cell   and  Strain  Gauge 
In  the  study  of  liquid  and  solid  He  a  measurement  of   pressure  can 
provide  very  fundamental   information  about  the  sample.    The  partition 
function  2  containing  the  Hamiltonian  of   the  microscopic  system  gives 
all  the  thermodynamic   quantities.    The  pressure  is  related  to  the 
partition  function  Z  by 

P  =  kBT   (3inZ/3V)T     .  (2.8) 

The   capacitance  strain  gauge,   proposed  first  by  Straty  and  Adams 
(1969),   measures   the  sample   pressure   in  situ.    A  sample  cell   and  a  strain 
gauge  used  in  this  experiment  are  designed  specifically  to  operate  at 
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1    mK  temperature  range  and  pressure  along  the  melting  curve  of  ^He  with 
a  sensitivity  AP  =  0.1    Pa  or  AP/P  =  3x10~8. 

The  -*He  sample  cell  contains  three  main  pieces  as  shown  in  Figure 
2.5.   The  part  B  is  a  flexible  diaphragm  which  under  the  pressure  change 
gives   the   different  capacitance  for  measurement.    The  part  A  and  B  form  a 
sample  chamber.    The  JHe  is  introduced  through  a  filling  capillary  into 
the  chamber    (not  shown  in  the  figure).    An  indium  "0"  ring  seals  the  gap 
between  A  and  B  very  well   at  the  melting  pressure  or  even  higher 
pressure  over   10  MPa. 

The  design  of   the  diaphragm  geometry  is  based  on  the  formula  given 
by  Straty  and  Adams   (1969) 


._     15Pm     3/2 


S     =   15a2P     /4t2    ,  (2.10) 

y  m 


where  a  is  the  diaphragm  radius,    t  the  diaphragm  thickness,    Pm  the 
maximum  pressure  in  the  measurement,    S     the  yield  stress  of  the 
material,    E  the  modulus  of  elasticity  of   the  material,    1   the  distance 
between  two  plates  of  the  capacitor,    AP  the  pressure  sensitivity 
required,    AC/C  the  relative  sensitivity  of   the  capacitance  bridge. 

The  sample  cell  sits  at  the  center  of  magnetic  field.   The  cell   is 
designed  to  meet  a  number  of  special  requirements  of  measuring  pressure 
as  a  function  of   temperature  and  magnetic  field  with  high  precision.    The 
cell   body  is   constructed  of  sterling  silver.    The   gauge  is  made  of 
beryllium  copper.    Eight   0.76  mm  diameter   pure  silver  wires  are  welded  to 
the  cell  body  to  provide  additional   contact  with  the   packed  silver 
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Figure  2.5       3He  sample  cell 

A:    heat  exchanger,      B:    diaphragm  and  movable  plate,      C:    fixed  plate 
holder,      1:   holes,      2:   Ag  wires  for  thermal   contact,      3:   Ag  packed 
powder,      4:    indium  seal,      5:    capacitor  plates,      6:   Pt  NMR  thermometer 
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powder.  Another  eight  1  mm  holes  drilled  along  the  length  of  the  silver 
powder  and  a  central  5  mm  diameter  hole  give  a  great  contact  between  the 
liquid  JHe  and  silver  powder.  The  powder  is  packed  from  3.44  g 
700  A  commercial  silver  manufactured  by  Vacuum  Metallurgical  Co.  LTD. 
(Japan)  at  25  MPa  and  left  unsintered.  The  filling  factor  after  the 
packing  is  40.6$.  The  surface  area  of  the  power  is  7.8  m2.  The  open 
space  of  the  cell  is  0.49  cm3  and  the  pore  space  of  the  packed  silver  is 
0.84  cm  .  All  these  means  have  greatly  improved  the  thermal  contact 
between  the  demagnetization  stage  and  helium  sample  in  the  cell. 
A  strain  gauge  located  out  of  the  magnetic  field  provides  a 
comparison  with  the  melting  pressure  in  the  sample  cell  sitting  at  the 
center  of  the  field.  The  gauge  is  made  of  beryllium-copper.  Only  0.61  g 
silver  powder  are  packed  in  the  gauge.  The  gauge  has  been  heat-treated 
at  350  C  for  30  minutes  to  improve  the  mechanical  stability  and  thermal 
property.  The  surface  area  of  the  sintered  powder  is  0.55  m2. 

2.7  Pt  Pulsed  NMR  Thermometer 
For  nuclear  paramagnets,  such  as  platinum  or  copper,  their  nuclear 
susceptibility  follows  Curie's  law 

C 
nuc 

Xnuc     T     '  (2.11) 

n 

for  the  temperature  higher  than  the  ordering  temperature,  where  y    is 

Anuc 

the  nuclear  susceptibility,  Cnuc  is  the  Curie  constant  for  the  nuclear 
spin  system,  and  Tn  is  the  nuclear  spin  temperature.  By  measuring  the 
nuclear  susceptibility  the  temperature  Tn  can  be  determined  through 
Equation  (2.11).  Because  the  nuclear  ordering  temperature  is  of 
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microkelvin  range,   this  sort  of  thermometry  is  very  suitable  for  any  low 
temperature  yet  achieved  in  bulk  matter. 

There  are  several  advantages  of  using  platinum  as   the  material  of 
the  thermometer.   Firstly,   platinum  has  relatively  long  spin-spin 
relaxation  time  t2   (1   ms)  which  simplifies   the  observation  and  recording 
of  the  NMR  signal.   Secondly,   the  spin-lattice  relaxation  time  t     follows 
Korringa's  relation  (1950) 

Ve  =  K     '  (2-12> 

where  Te  is  the  conduction  electron  temperature,  and  K  is  the  Korringa 
constant.  Because  the  Korringa  constant  of  platinum  is  very  small 
(0.0296  sec-K,  by  Aalto  et  al.,  1972),  the  T,  is  very  short  even  at  1  mK 
temperature  range.  This  results  in  rapid  equilibrium  between  the  nuclear 
spin  temperature  Tn  and  conduction  electron  temperature  T  .  Thirdly, 
there  is  only  one  magnetic  isotope  195Pt,  which  removes  the  beat 
structure  of  the  NMR  signal. 

The  pulsed  NMR  technique,  shown  in  Figure  2.6,  is  used  in  this  work 
to  measure  the  nuclear  susceptibility  of  195Pt.  In  a  steady  applied 
field  Bq  the  nuclear  magnetization  M  is  directed  parallel  to  BQ ;  M  may 
be  tipped  through  an  angle  9  by  applying  a  small  field  B1  perpendicrlar 
to  ^r,.  The  tipping  angle  depends  on  how  long  the  pulse  B1  lasts.  The 
pulse  could  be  chosen  as  a  finite  burst  of  sinusoidal  field  B^incot. 
After  the  tipping  pulse  the  nuclear  magnetization  precesses  around  the 
steady  field  1q  with  Larmor  frequency. 

w0  =  YB0  '  (2-13) 
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Figure  2.6        Principle  of  pulsed  NMR   thermometry 
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where  y  is  the  gyromagnetic  ratio  of  195Pt.   The  precessing  magnetization 
induces  in  a  receiver  coil  an  oscillating  signal  which  dies  away  as  the 
nuclear  spins  return  to  thermal  equilibrium.   This  free  induction  decay 
(FID)   is  characterised  by  time  x2    ,   the  spin-spin  relaxation  time. 

The  NMR  thermometer  must  be  constructed  with  finely  devided 
material   because  the  rf  field  only  penetrates  a  skin  depth  into  a  metal 
sample.   The  skin  depth  can  be  written  as 


1  1/2 

d   =  27(^  •  (2.14) 


where  p  is  the  electrical  conductivity  of  the  material,   u     is  the 
magnetic  permeability  in  vacuum,   and  v   is  the  rf  frequency.    The  brush  is 
made  of  platinum  wire  of  25  urn  diameter,   guaranteeing  the  sample  is 
totally  penetrating  by  rf  signal.    The  highest  working  frequency  of  our 
NMR  thermometer  is  4.5  MHz,   corresponding  to  a  skin  depth  of  45  wm  for 
platinum. 

The  procedure  for  making  Pt  brush  was  as  follows:    (1)   Use  200  Pt 
wires  of  25  urn  diameter   (99.999?  purity)   to  make  a  bundle.    (2)    Insert 
the  Pt  bundle  into  the  slit  of  the  platinum  block  holder.    (3)   Weld  the 
bundle  to  the  platinum  block  holder  in  an  argon  gas  enviroment.    (4) 
Insulate  the  wires   4-5  times  with  Acrylic  Spray  Coating  (Krylon  No. 
1303A)  on  the  Pt  wires  of  the  bundle.    (5)    Bundle  the  Pt  wires   tightly 
together  to  get  a  high  filling  factor.   Cut  the  bundle  end  to  give  a 
length  of  1    cm. 

The  rf  coil,   made  of  50   ym  diameter   copper  wire,   has   two  layers, 
each  of  which   is  60   turns.    The  total   length  of   the  coil   is  5  mm,   and  the 
diameter   of  the   coil   is  2  mm.    After  winding  the  first  layer  Stycast  1266 
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was  used  to  glue  it  down,   and  allowed  to  dry.   Then  the  second  layer  was 
wound,   and  the  Stycast  was  used  to  fix  the  position  of  the  coil.   The 
inductance  of  the  coil   is  6.4  pH   .    Finally,   the  Pt  brush  was  placed  into 
the  rf  coil,   and  a  slight  amount  of  GE  vanish  was  used  to  prevent 
vibration.   The  NMR  thermometer  was  screwed  to  the  silver  rod   (thermal 
link)  near  the  ^He  sample  cell. 

2.8     Superconducting  Magnets 

There  are  two  superconducting  magnets:    one  produces  the  magnetic 
field  for  the  nuclear  demagnetization,   the  other  produces  the  field  on 
the  •'He  sample  cell   and  Pt  NMR   thermometer. 

The  magnet  for  the  nuclear  demagnetization  is  a  home-made  one, 
built  by  Yokio  Morii    in  1982.    NbTi  superconducting  wire  "Supercon  279E- 
9B1A"  was   used  for  the  winding  of  the  solenoid.    The  main  solenoid  is 
15.24  cm  long,    5.4  cm  inner   diameter.    There  are  32   layers,    total   10416 
turns.    In  order  to  obtain  high  homogeneity  in  the  demagnetization 
region,   a  pair  of  compensation  coil  were  added  at  the  ends  of  the  main 
solenoid.    The  field-current  ratio  of  this  magnet  is  0.0790   T/Amp.    The 
maximum  field  is  5   Tesla. 

The  other  superconducting  magnet   providing  the  field  at  the  ^He 
cell  and  NMR  thermometer  is  a  commercial  one,   made   by  Nalorac  Cryogenic 
Corporation.    The  maximum  field  is  4.7   Tesla.    Four  superconducting  shims 
are  used  to   provide  high  field-homogeneity,   which  is  typically  0.3  ppm 
in  a  10  mm  diameter  and  15  mm  long  sample  region. 

The  main  field  for  the  nuclear  demagnetization  and  the  NMR  field 
should  be   put   in  opposite  so  that  they  will   be   in  the  same  direction  in 
the  reversed  demagnetization  region.    Decreasing  the  demagnetization 
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field  will  cause  decreasing  the  total  field  in  the  region  where  there  is 
lots  of  copper  and  produce  a  supplemental  demagnetization  effect. 

Superconducting  switches  were  installed  on  both  magnets  and  all  the 
shims  to  facilitate  extremely  low  drift  operation.  When  the  supercon- 
ducting magnet  is  working  under  persistant  mode,   the  trapped  field  is 
found  to  decay  with  an  initial   time  constant  of  1  O4  hours  due  to  the 
process  of  flux  creep   (File  and  Mills,    1963).   The  superconducting  switch 
can  be  thermally  actuated  by  controlling  the  current  in  a  heater  which 
is  close  to  the  switch. 

2.9     Electronics  for  Pressure  Measurement 
The  capacitance  of  the  strain  gauge  and  3He  cell   are  measured  by 
using  two  sets  of  capacitance  bridge.    Through  the  calibration  of 
capacitance  versus   pressure  performed  in  advance,    the  3He  melting 
pressure  could  be  determined  accurately. 

Figure  2.7  shows   the  schematic  diagram  of  the  electronics  for  the 
pressure  measurement.    A  lock-in  amplifier    (PAR5204)   serves   as  oscillator 
and  detector.    It  supplies  a  sine  wave  of  2A  kHz,    8  V       to  the  bridge 
through  an  isolation  transformer    (Model  Gertsch  ST-200C).    The  ratio  arms 
of  the  bridge  are  the  transformer  windings  of  AC  Ratio  Standard,   Model 
1011  A,    from  Eaton  Corporation.    The  AC  Ratio  Standard   consists  of  seven 
transformer  windings  and  seven  rotary  switches.    The  ratio  accuracy  of 
the  standard  is   based  upon  the   use  of  a  toroidal  autotransformer  which 
is  not  affected  by  age  or  enviromental   conditions.    The  resolution  of  the 
Ratio  Standard  is   1x10~7«    The  Cp  in  Figure  2.7  is   the  capacitance  of 
the     He  cell.    The  Cpef   is  a  fixed  reference   capacitor  which  is  made  of 
beryllium-copper  and  heat-treated  at  350  °C  in  an  H2  enviroment  for  2 
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hours.   The  capacitance  of  the  reference  capacitor  is  about  50  pF, 
comparable  to  that  of  the  capacitance  of  the  sample  cell.    The  output 
from  the  bridge  is  amplified  by  the  lock-in  and  recorded  continuously  on 
a  chart  recorder.   Also,   the  output  from  the  lock-in  is  recorded  by  the 
HP3421A  Data  Acquisition  Unit  and  stored  in  the  HP9845B  computer  along 
with  the  NMR  signal  every  10  to  30  minutes. 

The  other  capacitance  bridge  for  the  strain  gauge  out  of  the 
magnetic  field  is  basically  the  same  as   that  for  3He  cell,   except  that 

the  reference  voltage  used  is  a  sine  wave  of  3-5  kHz  and  5.4  V 

2-10 Electronics  for  Pt  Pulsed  NMR  Thermometer 

Figure  2.8  is  the  schematic  diagram  for  the  Pt   pulsed  NMR 
thermometer.    The  HP3325A  synthersizer  sends   a  continuous  sine  wave  to  an 
electronic  gate  controlled  by  a  burst   pulse  timer.    The  frequency  of  the 
sine  wave   is  set  at  the  Larmor  frequency  of   '95pt  at   tne  magnetic  field 
applied  on  the  thermometer  and  3He  cell.    The  burst  pulse  timer,   made  by 
the  electronic  shop  of  the  Physics  Department,    is  triggered  by  the  Sync, 
output   of  the  HP3325A.    The  burst   pulse   timer  has  two  functions.    One  is 
to  provide  a  pulse,   the  width  of  which  is 

1  =   nT     •  (2.15) 

where  n  is  the  number  of  cycles  of  the  sine  wave,   and  T   is  the  period. 
The  pulse  opens  and   closes   the  electronic   gate.    The  tipping  angle  9   of 
the  magnetization  M  could  be  adjusted  by  setting  different  number  of 
cycles  n.    The  other  function  of  the  burst   pulse  timer  is  to  trigger  the 
Nicolet  204-A  digital  oscilloscope  which  stores   the  free  induction  decay 
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of  the  magnetization.   The  burst  pulse  timer  can  work  in  two  modes, 
automatically  or  manually. 

The  signal  output  from  the  electronic  gate  is  a  sine  wave  of  n 
cycles.    The  signal   passes  through  a  RF  power  amplifier,   a  suitable 
attenuator  and  a  tuning  tank  to  the  rf  coil   inside  the  cryostat.    This  rf 
signal   produces   the  small  magnetic  field  B1    perpendicular  to  the  applied 
field  Bq   (see  Figure  2.6).   The  L,   C,   and  C2  compose  a  broad  band  tuning 
circuit.   The  peak  of  the  frequency  characteristic  curve  of  the  tuning 
circuit  should  be  adjusted  to  the  Larmor  frequency  corresponding  to  BQ. 
The  L  is  the  inductance  of  the  rf  coil,   which  is  the  transmitter   coil  as 
well  as   the  receiver  coil.    The  C-,    is  a  variable  capacitance  for 
adjusting  the  resonance  frequency  to  the  Larmor  frequency.    The  C2  is  the 
capacitance  of  the  coaxial  cable  in  the  cryostat.   The  highest  frequency 
for  this  NMR  thermometer  is  4.5  MHz,   which  is  limited  by  the  existance 
of  C2.    Two  silicon  diodes   in  the  tuning  tank  compose  a  gate  which  allows 
the  excitation  pulse  to  the  resonance  circuit  and  FID  signal  to  the 
preamplifier.    A  resistor  of  470   kfl  is  used  only  for  tuning  C,    to  obtain 
the  maximum  output  from  the  tuning  tank  at   the  Larmor  frequency. 

The  FID  signal,    amplified  by  a  low  noise  amplifier    (I6db)   and  a 
multi-stage  amplifier    (45db),    is  stored  in  the  memory  of  Nicolet  204-A. 
The  time  integral  of  the  FID  is  to  be  calculated  by  the  HP9845B,   and 
then  stored  in  a  data  file  for   analysis  later. 

2.11      Electronics  for  Superconducting  Magnets 
The  superconducting  magnet  for  the  nuclear  demagnetization  and  the 
NMR  magnet  share  a  main  power  supply  Kepco  ATE6-100M.    For  most  of  the 
time  in  the  experiment  the  main  power  supply  serves  as  the  current 
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source  for  the  nulear  demagnetization  magnet.   Only  when  the  magnetic 
field  on  the  JHe  cell  and  NMR  thermometer  needs  to  be  changed,   is  the 
main  power  supply  switched  to  the  NMR  magnet. 

Figure  2.9  shows  the  electronics  for  the  superconducting  magnets.    A 
resistor  of  0.01    n  and  a   voltmeter  are  used  to  monitor  the  current 
through  the  magnet.   The  magnetic  field  can  be  easily  calculated  by  the 
calibration  done  before. 

The  HP41-CV  calculator  serves  as  a  digital  data  bus,   which  produces 
digital  data  within  a  certain  time.    The  Kepco  programmer  SN488-1  22  is 
used  as   the  interface  between  HP41 -CV  and  ATE6-100M.    The  SN488-122 
converts  the  digital  data  from  the  bus  to  an  analog  voltage  output, 
which  serves  as   the  control   signal    input  for  ATE6-100M.    In  response  to 
the  control  signal   the  ATE6-100M  is   programmed  to  the  magnitude 
initially  commanded  by  the  data  bus   controller. 

Another  supplementary  power  supply  is  available  for  the   persistent 
heat  switch  of  the  NMR  magnet.    Anytime  when  the  NMR   field  needs   to  be 
changed,    the  heater   is  energized  by  about   150  mA   into  its   5  Q 
resistance.   The  superconducting  wire  which  is  close  to  the  heater 
changes  to  the  normal  state.    After   a  necessary  adjustment  for  the  field, 
the  heater   is  turned  off,   and  the  magnet   is  put   back  to  the  persistant 
mode. 

Appendix  C   gives   the  HP41-CV  program  for   producing  digital   data. 


>  o 

*  "3 

Q.  i» 


CHAPTER     3 
EXPERIMENTAL  TECHNIQUES 

3-1      Refrigeration 
The  techniques  for  refrigeration  used  in  the  experiment  are 
dilution  refrigeration  and  nuclear  demagnetization.    An  SHE  DRP-43 
dilution  refrigerator  runs   continuously  to  keep  the  cryostat  as  cold  as 
about  5  mK.   Unlike  the  dilution  refrigeration,   the  nuclear  demagneti- 
zation is  a  one  shot  operation.    The  period  of  one  cycle  lasts  about  2  to 
3  weeks  in  this  experiment. 

The  procedure  for   cooling  down  the  cryostat  from  room  temperature 
(300  K)   is  as  follows.    Before  the  precooling  by  liquid  N2,   all  the 
traps,    the     He  cell  and  strain  gauge  should  be   pumped  out  for   12  hours 
or  longer   at  room  temperature.    The  electrical   leads  for   all   the 
components,    including  the  capacitance  gauges,    NMR  coil,   superconducting 
magnets,   heat  switches,   carbon  resistors,   and  heaters  have  to  be  checked 
out  at  room  temperature.    Leak  checking  of  the  vacuum  can,    3He  cell, 
strain  gauge  and  circulation  system  of  the  dilution  refrigerator  must   be 
done  carefully.    If  everything  goes  normally,   then  it  is  ready  to  precool 
to  liquid  N2  temperature   (77  K).    The   vacuum  can  is  filled  with  dry  N2 
gas   to  1000  Mtorr  as   exchange  gas.    The  dewar  used  in  the  experiment   is 
superinsulated  and  has  no  liquid  N2  jacket.    The  cryostat   is  preceded 
with  liquid  N2  for  5   to  8  hours.    After  the  innermost   parts   of  the  dewar 
reach  77  K,    the  remaining  liquid  N2   should  be  transfered  out.    The 
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electrical   test  and  leak  test  are  repeated  at  the  liquid  N2  temperature. 
Exchange  gas     He  is  then  put  into  the  vacuum  can  at  a  pressure  of 
1000  utorr.    Liquid  4He  is  transfered  slowly  into  the  dewar  at  the 
beginning  in  order  to  use  ^He  efficiently.   When  the  temperature  reaches 
about  12  K,   the  transfer  process  should  be  interrupted  for  about  H  hours 
for  pumping  the  residue  Htie  exchange  gas  from  the  vacuum  can.   After  that 
the  transfer  is  resumed  until  the  dewar  is  full  of  liquid  4He. 

Under  normal  operating  conditions,   the  average  consumption  of 
liquid     He  is  22  liters   per  day.    Liquid  4He  is  transfered  every  three 
days  on  schedule. 

After  the  first   time  of   transfer,   liquid  HHe  in  the   1    K  pot   is 
pumped  by  Kinney  pump   (see  Figure  2.3).    The  temperature  of  the   1    K  pot 
maintains  at   about   1.3  K.    The   1    K  pot  fills  automatically  due  to  the 
siphon  effect. 

The  next  step  is  to  condense  the  3He/4He  mixture  into  refrigerator. 
At  first  the  gas  storage  tanks  are  opened  to  the  backside  of  the 
mechanical  pump  ED660.  When  a  part  of  mixture  gas  goes  into  circulation, 
and  the  pressure  of  the  tanks  drops  down  to  about  200  mtorr,  the  mixture 
gas  may  be  pumped  with  the  mechanical  pump  ED660  into  the  circulation 
until  all  the  gas  is  condensed  into  the  refrigerator.  This  process  takes 
about   4  to  6  hours. 

The     He  circulation  rate  is  controlled  by  a  heater   in  the  still. 
The  heater,   made  of  Karma  wire  with  0.076  mm  diameter  has  a  resistance 
of  533  ohms  at  room  temperature.   Usually  the  still  power   is  set  at   1    mW, 
giving  0.75  mmole/sec  ^He  circulation  rate.    This   power   could   be 
increased  to  2  mW  after  magnetizing  the  demagnetization  stage  in  order 
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to  speed  up  the  precooling  process.  It  takes  2  days  to  cool  the  mixing 
chamber  and  demagnetization  stage  down  to  7  mK  in  a  field  of  H   T. 

There  are  two  different  modes  for  the  nuclear  demagnetization.  The 
fast  mode  reduces  the  current  in  the  magnet  by  1  3  mA  per  step,  and  the 
slow  mode  reduces  the  current  by  1 . 3  mA  per  step.  The  time  interval 
between  steps  is  determined  by  the  operator.  The  nuclear  demagnetization 
is  typically  done  in  several  different  rate.  To  start  with  the 
demagnetization  a  relatively  fast  rate  is  chosen.  For  example,  0.1 
mT/sec  is  recommended.  After  the  field  is  reduced  from  4  T  down  to  about 
1.2  T,  the  rate  is  slowed  down  to  about  0.05  mT/sec.  With  this  rate  the 
liquid  He  in  the  cell  goes  through  the  superfluid  transition.  The 
demagnetization  is  hold  after  the  superfluid  transition  of  liquid  3He 
shows  up  on  the  chart  recorder.  It  takes  about  an  hour  to  reach  thermal 
equilibrium  for  the  3He  in  the  cell.  The  NMR  thermometer  is  calibrated 
in  the  vicinity  of  3He  superfluid  transition.  As  the  temperature  goes 
lower,  the  thermal  time  constant  gets  longer.  On  the  other  hand,  ramping 
the  field  down  would  cause  heating  due  to  eddy  currents  in  metal  in  the 
field  region.  Because  of  these  two  reasons  the  demagnetization  rate 
should  be  reduced  as  the  temperature  is  lowered.  Except  reducing  the 
demagnetization  rate  at  lower  temperature  there  are  several  pauses  of  2 
to  3  hours  for  each  in  the  demagnetization  process  to  ensure  the  melting 

o 

He  in  good  thermal  equilibrium  with  the  demagnetization  stage.  When  the 
temperature  is  close  to  the  ordering  transition  of  solid  3He,  the 
demagnetization  rate  is  controlled  at  5x10"^  mT/sec,  200  times  slower 
than  the  initial  rate.  After  slowly  cooling  through  the  ordering 
transition,  the  demagnetization  is  stopped.  The  cryostat  warms  up 
because  of  the  heat  leak  which  is  typically  1  nW.  The  warming  rate 
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varies  from  2.3  uK/hr  to  5.8  yK/hr,   depending  on  how  large  the  field 
left  on  the  demagnetization  stage.    The  oool   down-warm  up  process  usually 
is  repeated  twice  or  even  three  times  at  the  different  warming  rate  to 
see  the  effect  on  the  pressure  changing  as  a  function  of  time.    Each  cool 
down-warm  up  cycle  takes  a  few  days  to  a  week. 

When  the  cell   is  warmed  up  to  about   10  mK  which  corresponds   to 
almost  the  maximum  entropy  of  solid  3He,    Rln2   per  mole,    the  demagneti- 
zation stage  needs  to  be  magnetized  again. 

3-2     Pressure  Measurement 

The  calibration  of  capacitance  for  the  3He  cell  and  the  strain 
gauge  as  a  function  of  pressure  has   to  be  done  at   rather  high 
temperature  before  the  3He  in  the  cell  or  gauge  reaches  the  melting 
curve.    The  suitable  temperature  to  do  so  is  around  1    K. 

Figure  3.1    and  Figure  3.2  are  the  results  of  the  calibration  for 
the     He  cell  and  the  strain  gauge.    The  pressure  is  measured  by  a  HEISE 
pressure  gauge  which  has  an  accuracy  of  1x1 04   Pa.    The  capacitance  is  the 
reading  of   the  ratio  transformer,   rather  than  the  actual   capacitance 
value  of  the  cell  or  gauge.   We  fit  the  data  in  a  parabolic  function, 
having  a  form 


P  =   AQ   +  AlC  +  A2c2     ,  (3.1) 


where  AQ,    A1 ,    A2  are  the  fit   parameters,   and   c  is   the  capacitance. 
The  parameter  AQ   is  adjusted  for   fitting  the  A  transition  of 

■3 

superfluid  JHe.    The  capacitance  at  the  A  transition   could   be  determined 
very  accurately  on  the  chart.    Having  the  capacitance  corresponding  to 
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the  A  transition  from  the  chart  and  the  pressure  at  the  transition  based 
on  Greywall's  measurement   (1985),   the  constant  term  in  Equation   (3.1) 
can  be  readjusted  to  be  AQ     under  the  assumption  that  parameters  A1    and 
A2  remain  unchanged.   After  this  adjustment  a  new  calibration  of  pressure 
versus   capacitance  is  used  for  determining  the  melting  pressure, 


P  =   A0'+  A]C  +  A2c2      .  (3.2) 


In  a  magnetic  field  the  A   transition  of  superfluid  3He  splits   into 
two  separate  transitions   (Gully  et  al.,    1973).    For  example,   the 
separation  of   two  transitions  at  0.404   T  in  our  experiment   is  26  uK. 
This  splitting  increases  with  the  magnetic  field.    The  adjustment  for  the 
parameter  Aq     is  based  on  the  center   point   between  the  two  transitions. 

The  output  of   the  lock-in  amplifier  for  the  3He  cell   and  strain 
gauge  are  recorded  on  the  chart  recorder  continuously  and  by  the  data 
acquisition  unit  HP3421  A,    along  with  the  NMR   signal   at   certain  time 
intervals   (see  Figure  2.7).    The  output  of  the  lock-in  could  be  easily 
converted  into  the  capacitance  using  the  settings  of  the  lock-in  and  the 
chart  recorder.   Then,   by  using  the  Equation   (3.2)   the  melting  pressure 
can  be  determined. 

3.3  Thermometry 

3.3.1   Setting  the  Electronics 

The  Pt  pulsed  NMR  thermometer  used  in  this  experiment  is  attached 
on  the  silver  rod  near  the  3He  cell  located  at  the  center  of  the 
magnetic  field  provided  by  Nalorac  superconducting  magnet.  The  ratio  of 
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the  resonant  frequency  to  the  magnetic  field  for  Pt  metal   is  8.786 
MHz/T.   This  ratio  is  different  from  the  value  of  9.153  MHz/T  for  Pt 
nucleus  because  of  the  Knight  shift. 

Once  the  magnetic  field  on  the  3He  cell   is  decided,   the  frequency 
of  Pt  NMR  thermometer  is  also  determined  by  the  relation 

v  =  8.786   BQ    (MHz)      ,  (3.3) 

where  BQ   is  the  static  magnetic  field  on  the   3He  cell.    The  capacitor  C1 
in  the  tuning  tank  is  adjusted  so  that   the  peak  of  the  tuning  curve  is 
at  the  Pt  resonant  frequency  under  this  field.    The  tuning  procedure  is 
as  follows    (see  Figure  2.6).    A  continuous  sine  wave  from  HP3325A  is  sent 
directly  to  the  tuning  tank  through  a  resistor  of  470   kfl.    The  response 
to  the  different  frequency  describes  a  characteristic  tuning  curve,   such 
as  shown  in  Figure   3-3.    For   a  lower  resonant  frequency  corresponding  to 
a  lower  magnetic  field  the  tuning  curve   is  flatter  than  that  for  a 
higher  resonant  frequency.   This  is   because  a  larger   capacitance  C, 
causes  more  dissipation  and  reduces   the  Q  value  of   the  tuned  circuit. 
For  example,    the  Q  is  equal  to  19.4  for   a  field  of  0.430   T,    and  is   16.5 
for  a  field  of  0.266   T. 

The  amplitude  of  the  rf  signal   should  be   chosen  carefully  in  order 
to  avoid  saturation  of  the  NMR  signal   at  low  temperature.    For  the 
temperature  range  from  0.5  mK  to  10  mK,   which  this  experiment   is  most 
interested  in,    the   voltage  of  800  to   900  mV   (peak  to   peak)   proved 
suitable  for  the  rf  signal. 

It  is   important  to  realize  that  the  tipping  pulse  causes  heating  to 
the  nuclear  spin  system.    Since  in  the  Boltzmann  distribution  a  reduction 
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of  magnetization  from  M  to  Mcos8   is  equivalent  to  a  temperature  from  T 
to  T/cose,   the  temperature  change  for  a  small  tipping  angle  9   is 


AT=   c^sT  -TsITe2     '  (3'4) 


Because  the  heat  capacity  of  a  nuclear  dipole  system  C_   is  proportional 

p 

to  1  /T  ,  the  heat  introduced  by  the  tipping  pulse  is  inversely 

proportional  to  the  temperature  as  follows 


AT  C  a   02/T  .  (3.5) 


A  small   tipping  angle  9  should  be   used  in  order  to  reduce  this  heating. 
The  tipping  angle  is  determined  by  the  length  and  amplitude  of  the  rf 
pulse  as  follows 

9   -  2irn  =p     ,  (3.6) 

B0 

where  n   is  the  number  of  cycles,    B,    is  the  amplitude  of  rf  magnetic 
field,    Bq   is  the  applied  external  field.    In  our  experiment  n  was  75  and 
the  peak  to  peak  voltage  was  about  880  mV  which  corresponds  to  4x10~^T 
for  B1 .    For  the  external  field  BQ  from  0.266  T  to  0.495   T  covered  in 
this  work,    the  tipping  angle  is  from  2.7°  to  5°. 

3.3.2     Calibration 

The  NMR  thermometer   is  calibrated  against  the  superfluid  transition 
of  liquid  ■'He.    According  to  Greywall's  temperature  scale   (1985)   this 
transition  happens  at  2.708  mK  and  melting  pressure  3M.338  bar.   Appendix 
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D  is  the  program  for  calculating  the  temperature  from  melting  pressure 
based  on  Greywall's  scale  (Greywall,  1985). 

In  fact,  the  calibration  is  performed  in  the  vicinity  of  the 
superfluid  transition  of  ^He.  By  measuring  the  temperature  from  the 
melting  pressure  and  the  susceptibility  from  the  NMR  signal,  the  Curie 
constant  could  be  determined  by 

C  =  XT  •  (3-7) 

This  calibration  is  repeated  at  several  points  near  the  A  transition 
point  and  the  average  Curie  constant  is  taken  for  final  use. 

Because  the  field  for  the  nuclear  demagnetization  and  the  NMR  field 
are  not  totally  isolated  from  each  other,  the  field  change  at  the  demag- 
netization stage  slightly  shifts  the  NMR  frequency.  The  ratio  of  the  NMR 
frequency  shift  to  the  field  change  at  the  demagnetization  stage  is 
about  -10  kHz/T.  The  total  shift  of  the  NMR  frequency  in  the  whole  de- 
magnetization process  from  4  T  to  almost  zero  field  is  40  kHz,  appro- 
ximately 1%   of  the  NMR  frequency.  Whenever  the  field  at  the  demagneti- 
zation stage  is  changed,  the  rf  frequency  needs  to  be  adjusted  to  match 
the  Larmor  frequency.  The  NMR  frequency  could  be  adjusted  by  using 
another  rf  signal  to  beat  against  the  NMR  signal  or  by  sweeping  the  rf 
frequency  to  find  the  maximum  NMR  signal. 

Data  are  taken  on  both  cool  down  and  warm  up.  But  the  data  from  the 
warm  up  is  more  reliable  because  the  field  at  the  demagnetization  stage 
and  NMR  frequency  are  not  changing  in  the  warm  up  process.  The  Curie 
constant  obtained  at  warm  up  is  used  for  calculating  the  temperature. 
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3.H     Data  Acquisition 
The  information  provided  in  the  experiment  is  recorded  continuously 
by  a  three-pen  chart  recorder  and  by  a  HP9845B  computer  digitally  point 
by  point. 

The  chart  recorder  describes  the  pressure  change  in  the  3He  cell 
and  the  strain  gauge  as  a  function  of  time.    In  the  process  of  nuclear 
demagnetization  the  chart  recorder  also  records  the  current  in  the 
magnet.   The  HP9845B  computer   collects  the  data  from  HP3421A  data  acqui- 
sition unit  and  the  Nicolet  digital  oscilloscope  204-A.    There  are  two 
modes  for  the  data  taking,   automatic  and  manual.    In  the  automatic  mode 
the  burst  pulse  timer  sends  a  pulse  to  trigger  the  NMR  circuit  at 
intervals  which  are  adjustable  from  a  few  minites  to  an  hour.    This  mode 
is  useful  for  monitoring  the  temperature  and  pressure  change  over  a  long 
time  period  of  a  couple  of  days  in  a  cool  down  or  warm  up  process.    In 
the  manual  mode  the  burst  pulse  timer  sends  a  pulse  only  at  the 
operator's  request.    This  mode  is  often   used  in  calibrating  the  NMR 
thermometer.    The  FID  signal    is  stored  in  the  memory  of  Nicolet  204-A. 
After  the  NMR  circuit  is  triggerd,    the  HP9845B  will   command  the  Nicolet 
204-A  to  transfer  the  data  of   the  FID  signal   to  the   computer  and 
calculate  the  time  integral  of  the  FID.    Also  the  HP9845B  asks  the 
HP3421A  data  acquisition  unit  to  scan  the  lock-in  outputs  which  could  be 
converted  to  the  capacitance  of  the  3He  cell  or  the  strain  gauge 
according  to  the  lock-in  settings.    All  the   information  about  the  Pt 
susceptibility,   lock-in  outputs,   and  the  time  provided  by  the   clock 
HP98035A  is  permanently  stored  onto  a  magnetic  tape  for  future  analysis 
and  printed  out  as  a  reference. 
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The  data  acquisition  is  controlled  by  a  program  "TEMP5A",  which  is 
shown  in  Appendix  E. 


CHAPTER     4 
DATA   REDUCTION 

Data  have  been  taken  in  magnetic  fields  of  0.266,    0.373,    0.390, 
0.400,    0.402,    0.404,    0.410  and  0.495  T  over  a  period  of  10  monthes.    For 
the  fields  of  0.266,    0.373  and  0.495   T  solid  3He  has   been  cooled  down  to 
about  0.5  mK,    and  then  warmed  up  to  10  or  20  mK.    Each  cycle  of   cool  down 
and  warm  up  took  about  one  month.    For  the  other  fields  we  only  concen- 
trated on  the  region  of  the  ordering  transition. 

4.1      Temperature  Smoothing 
Because  the  scatter   in  Pt  susceptibility  measurement  is  around 
0.3?,   the  individual  points  are  not  used  for  determining  the 
temperature.    Instead,    the  warming  rate  within  a  certain  time  is  first 
found  by  fitting  the  temperature-time  relation.    This  relation  is  linear 
in  a  short  time  period  of  one  day  or  even  longer.    But,    for   the  warm  up 
covering  2  or   3  days   the   best  fit  would   be   parabolic.    The  temperature 
assigned  to  each  datum  is  determined  by  the  fit  function  T(t)   based  on 
time  t,   which  is  recorded  on  a  magnetic  tape  along  with  the  other 
information  by  the  HP9845B  computer.    This  procedure  gives  smooth 
temperatures  required  for  taking  the  derivative  dP/dT  from  P(T).    Figure 
4.1a  is  a  linear  fit  in  a  40  hour   period  for  the  field  0.266   T.    The 
warming  rate  is  approximate  2.9   yK/hr.    Figure  4.1b  shows  a   better  fit  by 
using  a  parabolic  function  for  the  same  field  and  time  interval  as 
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Figure  4.1a.    For  calculating  the  temperature  within  a  long  time  period, 
a  parabolic  fit  between  time  and  temperature  should  be  used.    For  a  short 
time  period  about   1    day  a  linear  fit  is  usually  used  to  calculate  the 
temperature. 

M.2     Pressure  as  a  Function  of  Time  or  Temperature 
To  obtain  the  pressure  as  a  function  of  time  or  temperature  is 
quite  straight  forward.   A  data  file  contains  the  lock-in  outputs  and  the 
time  when  the  data  are  taken.   The  capacitance  of  the  %e  cell  or  strain 
gauge  could  be  converted  from  the  lock-in  output  by  using  the 
sensitivity  AC/AV  which  is  the  ratio  of  capacitance  change  relative  to 
the  lock-in  output.   Then  the  pressure  corresponding  to  that  capacitance 
is  calculated  by  Equation   (3-2). 

In  a  long  time  period  the  ratio  transformer  has   to  be  adjusted  in 
order  to  trace  all  the  pressure  change  on  the  chart  recorder.    It  happens 
sometimes   that  this  adjustment  of   the  ratio  transformer   could  cause  a 
small   discontinuity  on  the  pressure-time  curve.    A  constant  shift  for  the 
piece  which  does  not  match  the  whole  curve  usually  gives  a  satisfactory 
smooth  curve  of  pressure  versus   time.    This   procedure  is   also  necessary 
for  taking  the  derivative  dP/dT  from  P(T). 

From  the  pressure  as  a  function  of  time  P(t),   obtained  from  the 
above  procedure  and  the  fit  of  T(t)    it   is  easy  to  convert  the  P(t)   into 
P(T),    the   pressure  as   a  function  of  temperature.    Bacause  of  the  Kapitza 
resistance  between  the     He  sample  and  the  cell   body,    there  is  a 
temperature  discrepancy  between  ^He  and  the   cell   body   (or  the  nuclear 
demagnetization  stage).    The  thermal   time  constant   is  about   1    hour   around 
1    mK.    The  warming  rate  is   controlled  at  2  to  6   uK/hr.    The  temperature  of 
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He  thus  lags  about  2  to  6  yK  behind  the  cell  body.   On  the  other  hand 
this  temperature  gradient  will  cause  a  small  error  for  taking  dP/dT  in 
the  vicinity  of  the  ordering  temperaure,   which  is  of  little  consequence 
here  since  we  make  no  use  of  the  temperature  dependence  of  solid  entropy 
derived  from  dP/dT. 

After  the  ordering  transition  happens,   the  temperature  of  the 
nuclear  demagnetization  stage  and  3He  cell  body  keep  the  same  warming 
rate  as  before  because  the  heat  capacity  of  the  nuclear  demagnetization 
stage  is   very  large  and  the  latent  heat  absorbed  by  solid  3He  at  the 
transition  is  very  small. 

Figure  4.2  shows   the  pressure  versus   cell   temperature  at  0.373   T. 
The   pressure  is  measured  relative  to  the   pressure  of  A  transition  P.. 
Figure  4.2a  covers  a  wide  temperature  range   up  to  10  mK.    The  detail 
around  the  ordering  transition  is  shown  in  Figure  4.2b,    where  the  dots 
are  the  data  points  and  the  dashed  line  illustrates   the  behavior  of  P(T) 
during  the  transition. 

The  results  for  other  fields  will  be  shown  in  the  next  chapter. 

4.3     Entropy  of  Solid  3He 
Based  on  the  pressure  as  a  function  of  temperature  obtained  from 
the  previous   analysis  the  entropy  of  solid  3He  can  be   calculated  by 
applying  the  Clausius-Clapeyron  equation 


dP       Ss   -Sl 

dT        V      -  V  (1  .31  ) 

s         1 


From  Equation   (1 .31  ) 
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sQ  =  J?  (v  -  V  )   +  S.     ,  (H.1) 

s        dT        s       1  1 


where   (Vg-V-^)   is  -1.31    cm^/mole  at  temperature  range  0-20  mK   (Scribner 
et  al.,    1969;  Halperin  et  al.,    1978).    Liquid  ^He  has  very  small  entropy 
in  its  Fermi   liquid  regime   (Wheatley,    1975;   Varma  et  al.,    1976) 


S,    =   YRT  for     T   «  T„      ,  (4.2) 

1  r 


where  Y  =  4.6  K   near  the  melting  volume,  Tp  is  the  Fermi  temperature, 

which  is  around  1  K  depending  on  the  pressure  (Leggett,  1975;  Wheatley, 

1975).  At  about  2.7  mK  liquid  -*He  undergoes  a  superfluid  transition  and 

has  an  even  lower  entropy.  For  the  purpose  of  calculating  the  entropy  of 
solid  JHe,  the  liquid  entropy  below  the  superfluid  transition  can  be 
neglected. 

The  derivative  dP/dT  is  taken  by 


dP  (J   +  M/91     P(N  +  I)  -  P(I)  ,,  _. 

-  (I  +  N/2)  =  T(N  -  I}  _  T(I)    ,  (4.3) 


where  P(I),    T(I)   refer  to  the  pressure  and  temperature  of  the  Ith  datum 
point;   N  is  an  even  integer.    A  suitable  N  is  chosen  based  on  the  density 
of  the  poits.   Usually  N   is  6  to  10,    corresponding  to  a  time  interval  of 
about  an  hour.    Before  taking  the  derivative,    the  pressure  values  are 
smoothed  by  a  subroutine  in  a  program  "SPLOT"   (see  Appendix  F) .    The 
value  of  dP/dT  by   Equation   (4.3)    corresponds  to  the  temperature  of 
[T(N+I)    +  T(I)]/2. 
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We  have  tried  another  way  to  find  the  derivative  dP/dT.  First,  fit 
N  points  by  a  polynomial,  then  take  the  derivative  analytically.  This 
method  does  not  improve  the  scatter  of  the  entropy  for  solid  ^He 
significantly. 

The  method  based  on  Equation  (^.3)  has  been  used  for  calculating 
the  entropy  of  solid  ^He  at  differen  magnetic  fields.  But  it  is  not 
applicable  for  the  data  points  around  the  ordering  transition  because 
the  temperature  of  ^He  is  not  in  thermal  equilibrium  with  the  cell  body. 

Appendix  F  gives  the  program  for  calculating  dP/dT  and  plotting  the 
entropy  S(T) . 


CHAPTER     5 
RESULTS  AND  DISCUSSION 

U — General  Features  for  Different  Order  Phase  Transition 
A  first-order  phase  transition  is  defined  by  a  discontinuity  of  the 
first-order  derivative  of  the  Gibbs  function  at   the  phase  transition. 
For  a  second-order   phase  transition  the  first-order   derivatives  of   the 
Gibbs  function  remain  unchanged,   but  the  second-order  derivatives  of  the 
Gibbs  function  undergo  finite  changes.    Figure  5.1   shows  the  main 
features  of  the  Gibbs  function,   entropy   (the  first-order   derivative  of 
Gibbs  function  respect  to  temperature),   and  specific  heat    (proportional 
to  the  second  derivative  of  Gibbs  function)   for  the  first-order  and  the 
second-order   phase  transition.    The  superfluid  transition  of  liquid  3He 
is  an  example  of  first-order  transition.    An  example  of  second-order 
transition  is  the  transition  from  superconductor  state  to  normal  state 
in   zero  magnetic  field. 

The  most  interesting  higher-order   phase  transition  is  lambda 
transition  which  is  accounted  for  by  the  fact  that  the  shape  of  specific 
heat   versus   temperature  curve  resembles   the  Greek  letter  lambda   (Figure 
5.1c).    The  transition  from  ordinary  liquid  UHe   (He   I)   to  superfluid 
liquid   (He  II)    is  a  typical   example  among  many  lambda  transitions. 
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(a) 


(b) 


(c) 


Figure  5.1        Features   of  Gibbs  function   (G),   entropy    (S),   and 

specific  heat    (C)   for   different  order   phase  transition 

(a)  first-order  transition 

(b)  second-order  transition 

(c)  X  transition 
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5.2     Ordering  Transition  of  Solid  3He  at  Melting  Pressure 
It  i3  important  to  understand  the  dynamical  behavior  of  P(t)  as 
distinct  from  P(T)   in  order  to  indicate  the  order  of  the  various 
transitions  in  this  experiment.   Because  of  the  very  large  heat  capacity 
of  the  nuclear  demagnetization  stage  and  small  latent  heat  absorbed  by 
solid  JHe,   the  warming  rate  of   the  cell   body  is  unaffected  by  the 
entropy  taken  up  by  solid  3He  in  the  phase  transition.   The  pressure 
measured  by  the  3He  cell  reflects  the  temperature  at  the  liquid-solid 
interface.    Thus  during  a  first-order  transition  there  is   a  sloping 
platean  in  P(t),   as   the  interface  warms  up  while  the  interior  of  the 
solid  undergoes  the  transition  (Figure  5.2a).    The  shape  of  the  plateau 
region  depends  on  the  warming  rate.    A  slow  warm  up  has  a  flatter  plateau 
than  that  for  a  fast  warm  up.    A  second-order  transition  with  a  specific 
heat  discontinuity  will  produce  a  discontinuity  in  dP/dt  although  dP/dT 
is   continuous,    implying  that  the  entropy  of  solid  3He  is   continuous 
(Figure  5.2b).    A  rapid  change   in  dP/dt  will   occur   in  the  case  of 
a  A  type  transition   (Figure  5.2c).    This  "enhancement"   in  the  behavior  of 
dP/dt  relative  to  dP/dT   is  well   known  from  the  observation  of  the 
superfluid  transition  in  liquid  3He. 

Without  analyzing  the  data  the  basic  feature  of  solid  3He  ordering 
transition  could  be  seen  clearly  on  the  chart. 

Figure  4.2  shows   the  pressure  versus   cell   temperature  at  a  field  of 
0.373   T.    Figure  4.2a   covers  the  temperature  range  up  to  10  mK.    Figure 
4.2b  shows  the  detail   around  the  transition.    The  plateau  region  reflects 
the  nature  of  first-order  transition.    The  cell   temperature  is   propor- 
tional  to  the  time  scale.    The   point  A  marks   the  onset  of  the  first-order 
transition.    Between  points  A  and  B  a  temperature   gradient  AT  =   10  yK 
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(a)        1st    order 


time 


(b)        2nd   order 


time 


(c) 


time 


Figure  5.2       P(T)   and  dynamic  behavior  of  P(t)   for 
different  order   phase  transition 

(a)  first-order  transition 

(b)  second-order  transition 

(c)  \  transition 
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accumulates  due  to  the  latent  heat  absorbed  by  the  solid  ^He  in  the 
phase  transition.   During  this  interval  and  a  few  hours  afterwards  the 
JHe  and  thermometer  are  far  from  equilibrium.    The  behavior  of  P(T)    in 
the  vicinity  of  the  transition  is  illustrated  schematically  by  a  dashed 
line.   Elsewhere  we  convert  P(t)   to  P(T)   by  using  the  measured  warming 
rate  dT/dt. 

The  other  runs  at  magnetic  fields  of  0.266  T  and  0.390  T  have  a 
similar   plateau  region.    Figure  5.3a  is  the  pressure  measurement   covering 
up  to  10  mK  for   0.266   T.    Figure  5.3b  and  Figure  5.4  are  the  details  of 
the  transition  region  for  0.266  T  and  0.390   T. 

Figure  5.5  shows   the  pressure  versus   cell   temperature   (or  time)  at 
a  field  of  0.495   T.    The  detail  around  the  transition  is  shown  in  Figure 
5.5b  where  the  point  C  marks  the  transition  point.   A  plateau  region  of 
P(t)  has  not  been  seen  within  the  resolution  of  this  experiment, 
implying  that  the  transition  is  not  first  order.   The  another  run  at  the 
magnetic  field  of  0.465  T  shows  similar   behavior  to  that  of  0.495  T 
(Figure  5.6).    For  the  HFP-PP  transition,   dP/dt  always  show  the  rapid 
change  characteristic  of  a  A  transition.    This  is   particularly  apparent 
at  fields  not  too  near  0.400  T,   so  that  the  transition  region  is 
broader.   A  similar  appearence  would  be  also  caused  by  a  second-order 
transition  with  a  temperature  gradient.  We  conclude  that  the  HFP-PP 
transition  is  not  first  order  for  0.4   <  B  <  0.495   T.    However,   we  cannot 
exclude  a  transition  with  a  small  entropy  discontinuity,   AS/Rln2  £  0.05, 
which  would  produce  an  indiscenrnible  plateau  in  P(t). 

In  the  warm  up  at  the  magnetic  field  of  0.400  T  there  are  two 
transitions  which  display  different  features    (Figure  5.7).    The 
transition  at  lower  temperature  shows  a  plateau  although  this  plateau 
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lasts  for  a  relatively  short  time  compared  with  the  plateau  at  low 
field.   The  other  transition  happening  at  higher  temperature  has  a 
similar   behavior  to  that  at  higher  field  of  0.495   T.    This  indicates  that 
the  transition  at  lower  temperature  is   a  first-order  transition,    the 
other  one  at  higher  temperature  is  not  a  first-order  transition.    Several 
other  runs   at  the  fields   0.402,    0.404,    0.410   T  also  have  the  same 
feature  as  0.400  T.    Figure  5.8  through  Figure  5.10  show  P(T)    for  these 
three  fields. 

To  summerize  what  we  have  seen,   we  understand 

1.  The  transition  is  first  order  only  when  one  of  the  phases  is 
LFP,    for  example,   the  transition  between  the  LFP  and  PP  at  low  fields  or 
the  transition  between  LFP  and  HFP. 

2.  The  HFP-PP  transition  is  not  first  order  for  0.4  <  B   <  0.495   T. 
It  is   possibly  a  second-order  transition  or  a  X  transition. 

Upon  taking  the  derivative  dP/dT  and  substituting  all  the  other 
values  for  the  quantities  in  the  Clausius-Clapeyron  equation  as 
discussed  in  section  4.3,   we  have  found  the  entropy  of  solid  ^He  as  a 
function  of  temperature.    Figure  5.11    shows  the  molar   entropy  of  solid 
He  for  0.266   T.    The  entropy  increases  abruptly  at  the  ordering 
temperature,   and  then  increases  gradually,   finally  reaches   its  maximum 
value  Rln2  at  high  temperature  around  10  mK.    The  insert  shows  the  detail 
of  the  entropy  change  near  the  transition  where  there  is  an  entropy 
discontinuity  ASg  =  0.35  Rln2.    Figure  5.12  gives  the  result  for  0.373   T. 
It  shows  a  similar  feature  as  Figure  5.11    for  0.266  T,    except  the 
ordering  transition  happens  at  a  little  lower  temperature,   and  the 
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entropy  discontinuity  is  about  AS     =  0.26  Rln2.   When  the  field  goes 
higher  this  entropy  discontinuity  at  first-order  transition  gets 
smaller.   At  0.390   T,    ASg  =  0.23  Rln2.   Osheroff  and  Yu  (1980)  reported 
that  ASg  =  0.44   Rln2  for  the  field  0.014  T.    Our  result  and  Osheroff  s 
result  are  qualitatively  consistent.   Table  3  lists  the  ordering 
temperature  and  the  entropy  discontinuity  of  the  first-order  transition 
below  0.400  T  from  our  measurement. 


Table     3 

Ordering  temperature  and  entropy  discontinuity  of 
the  first-order  transition  below  0.390   T 


B   (T)  Ordering  temperature   (mK)  AS  /Rln2 

s 

0.014*  1.028  0.443 

0.266  1.003  0.35 

0.373  0.943  0.26 

0.390  0.895  0.23 


It  is   difficult  to  make  an  accurate  entropy  analysis  for  the  case 
which  covers   two  transitions  from  0.400  T  to  0.410  T.    Table  4  lists  the 
transition  temperatures  and  the  interval   between  these  two  transitions 
for  the  field  from  0.400  T  to  0.410  T. 


*  From  Osheroff  and  Yu   (1980) 
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Table     4 
Transition  temperatures  for  fields  0.400  T  -  0.410   T 

B  (T)  Tc1    (mK)  Tc2    (raK)  AT  (yK) 


0.400  0.875  0.905  30 

0.402  0.885  0.930  45 

0.404  0.880  0.945  65 

0.410  0.830  0.930  100 


The  interval  between  two  transitions  is  from  30  to  100  pK.    It  takes 
about   10  hours  to  one  day  to  go  through  two  transitions  if  the  warming 
rate  is  about  3   yK/hr.   Because  the  -*He  in  the  cell  and  the  thermometer 
are  not  in  good  thermal  equilibrium  when  the  first-order  transition 
happens  and  few  hours  later,   the  method  described  in  section  4.3  for 
calculating  the  entropy  is  not  suitable  for  this  situation. 
Nevertheless,   the  nature  of  these  two  transitions  is  quite  clear  from 
the  trace  of  pressure  as  a  function  of  time  or  cell  temperature.   An 
approximate  value  of  the  entropy  discontinuity  of  the  first-order 
transition  at  0.400  T  is  0.12  Rln2.    This  small   entropy  decrease  is 
because  most  of  the  entropy  has  already  removed  by  ordering  in  the  high- 
field  phase   (see  next  section). 

There  is  no  entropy  discontinuity  for  the  fields  of  0.465   T  and 
0.495   T.    Figure  5.13  shows  the  entropy  versus   cell  temperature  at 
0.495   T,    the   insert   is  the  detail  of   the  entropy  change.    The  ordering 
temperature  T     =  1.03  mK.    The  entropy  of  solid     He  increases   gradually 
before  the  temperature  reaches  TQ.    Once  the  temperature  passes   through 
T   ,    the  slope  of  the  entropy  curve  decreases   greatly.    If  we  take  the 
derivative  dS3/dT,   we  can  find  the  specific  heat   for  solid  ^He  by  using 
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Figure  5.11  shows  the  result  of  Cy  in  unit  Rln2  for  field  0.495   T. 

5.3     Magnetic  Phase  Diagram 
Figure  5.15  is  a  modified  magnetic  phase  diagram  of  solid  3He  based 
on  the  previors  phase  diagram  Figure  1.3  and  the  results  of  this  work. 
Below  0.390   T  we  have  only  seen  one  first-order  transition,   correspon- 
ding to  one  phase  boundary  between  the  low-field  phase  and  the 
paramagnetic  phase.    From  0.400  T  to  0.410  T  there  exist  two  different 
transitions,   corresponding  to  two  phase  boundaries  in  this  region,  one 
of  which  is  the  phase  boundary  between  the  low-field  phase  and  the  high- 
field  phase,   and  the  other  of  which  is  the  phase   boundary  between  the 
high-field  phase  and  the  paramagnetic  phase.    For  the  higher  fields   the 
transition  between  the  high-field  phase  and  the  paramagnetic  phase  is 
not  first  order.   The  transition  temperature  increases  with  the  increase 
of  magnetic  field. 

The  open  circles   in  Figure  5.15  are  obtained  by  warming  up  the 
sample  through  the  transition  at   the  fixed  magnetic  fields.    Not   all  the 
results  are  shown  in  the  figure  in  order  to  make  the   intersection  region 
more   clear.    The  solid   circles  are  the  results   of  Osheroff ,    Cross  and 
Fisher    (1980).    The  dashed  lines  are  from  the  measurement   by  Adams  et  al. 
(1980).    The  solid  line  is   the  fit   partially  based  on  our  data  and 
partially  based  on  Osheroff's  data  below  0.8  mK.    This  solid  line  defines 
the  phase  boundary  between  the  low-field  phase  and  the  paramagnetic 
phase   (LFP-PP),   and  the  phase  boundary  between  the  low-field  phase  and 
the  high-field  phase   (LFP-HFP).    The  phase  transition  across  this 
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Figure  5.15       Magnetic  phase  diagram  of  solid  3He   (modified) 

Open  circles  are  from  this  work.    Solid   circles,    from  Osheroff 
Cross  and  Fisher    (1980).    Dashed  line,    from  Adams  et.   al.(1980).    Solid 
line  is  first-order  transition  line.    Dot-dashed  line  is  not  first-order 
line. 
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boundary  is  of  first  order.   The  dot-dashed  line  is  the  phase  boundary 
between  HFP  and  PP.   The  transition  across  this  boundary  is  likely  a 
second-order  transition,   or  perhaps  a  X  transition,   rather  than  a  first- 
order  transition. 

The  intersection  point  of  the  two  phase  boundaries  can  be  found  by 
the  following  way.    For  the  magnetic  fields  from  0.400  T  to  0.410  T  two 
different  type  phase  transitions  have  been  seen.   The  separation  between 
the  two  transitions  is  increasing  with  field  increase   (see  Table  4). 
Assuming  a  simple  linear  relation  between  the  separation  AT  and  the 
magnetic  field  B,   the  magnetic  field  at  the  intersection  point 
where  AT  =  0   is  determined  to  be  0.396   T  by  a  fit   (Figure  5.16).   The 
temperature  corresponding  to  this  field  is  0.90  mK. 

Compared  with  the  previors  phase  diagram  Figure  1.3  where  dB/dT  has 
a  discontinuity  along  the  LPF  boundary,   dB/dT   in  the  new  phase  diagram 
Figure  5.16  is  shown  as  continuous  along  the  LPF  boudary,    but  with  a 
rapid  change  near  the  intersection  point  where  the  PP-HFP   transition 
line  meets  the  first-order  transition  line.   This  is  similar  to  the 
pressure-temperature  diagram  of     He  where  the  A   line  meets   the  melting 


5.4     Spin  Waves  for   the  High-Field  Phase 
The  low  temperature  form  of   the  specific  heat  of   solid  ^He  gives 
important  information  on  the  magnetic  structure  of  the  high-field  phase 
or,    indeed,    if   it  is  simply  a  highly  polarized  paramagnet. 

For  the  temperature  range  T  €  Tc  the  temperature  dependence  of  the 
specific  heat  is  determined  by  the  low  energy  elementary  excitations  of 
the  system.    The   direct  interaction  between  nuclear  magnetic  dipole  could 
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be  ignored  at  millikelvin  temperature.  The  only  important  elementary 
excitations  are  the  spin  waves.  The  characteristic  magnitude  of  the 
spin-wave  frequency  is  kBTc/ft  and  the  spin-wave  velocity  C  is  of  order 
kBTca/fi  '  wnere  a  i3  the  lattice  constant. 

The  low  temperature  specific  heat  per  unit  volume  is  related  to  the 
combination  of  the  spin-wave  velocities  (Cross  and  Fisher,  1985) 

.     Ku  (k)  2   Ku  (k) 

where  n(x)  =  1/(ex-1)  is  Bose  distribution  function,  the  prime  denotes 
the  differentiation,  the  sum  runs  over  wavevectors  in  the  first 
Brillouin  zone  and  number  of  different  modes.  At  low  temperatures  the 
spin-wave  frequency  is  linear  in  the  wavevector  k.  By  using  the  linear 
expansion  for  small  k,  Equation  (5.2)  becomes 

o  2  i  4  »3 
17  kB  T       1 

CV  "  "J"  I  <— 3 >  .  (5.3) 

15  n°   a   CJ(k) 
a 

where  the  average  is  taken  over  angles  which  for  the  tetragonal  case  is 
given  by 

•   1   N  1 

<~^>  =  — 5 •  (5.4) 

C  J    C  1  C  l 
a      aj_  a  I 

By  fitting  the  pressure  below  the  TQ  to  the  T4  Osheroff  and  Yu  (1980) 
found  that  the  spin-wave  velocity  is  8.4  cm/sec,  which  is  roughly  of  the 
expected  magnitude. 

It  is  known  that  the  high-field  phase  has  a  high  susceptibility.  Is 
this  phase  a  magnetically  ordered  phase  or  just  a  highly  polarized 
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paramagnet?  The  partition  function  of  a  paramagnetic  system  i: 


m=+J      gyRBm 
Z  =  I   exp  C-j— y~)   ,  (5.5) 

m=-J      KB 


In  Z  =  in  sinh  (J  +  -)  a  -  In  sinh  |  a     ,  (5.6) 

where  J   is  the  spin,   g  Lande  g  factor,   yB  Bohr  nuclear  magneton,   B 
magnetic  field  and  a  =  gUgB/kgT.    If  the  solid  3He  at   high-field  phase 
were  really  a  paramagnetic  phase,   the  pressure  would  be 

P  =  PQ   +  AT  In   [exp   (0.775   |)    +  exp    (-0.775   |) ]      ,  (5.7) 

where  B   is  in  Tesla,    T   in  mK   (see  Appendix  B  for  derivation). 

Figure  5.17  shows   the  fit  to  the   data  points    (solid  line)   and  the 
pressure  of  a  free  spin  system   (dashed  line)    for  the  magnetic  field 
0.495    T.    The  dashed  line  is  a  very  poor  fit  to  the  data  points  because 
the  specific  heat  maximum  for  the  Schottky  anomaly  is   at   T  =  0.33  mK  and 
the  pressure  covered  in  the  temperature  range  0.5-1.0  mK   is  almost  into 
the  linear  region.    If  we  take  an  effective  field  as   a  fit   parameter  in 
P(T)    for   free  spins,   we  find  Beff  =   1.96   T   (Figure  5.18)   with  a   fit 
almost  as  good  as  T4  as  shown  in  Figure   5.17.    But  from  the  NMR 
experiment  reported  by  Adams  et  al.    (1980)    and  Osheroff  et  al.    (1980) 
the  NMR  frequency  shift  in  high-field  phase  is   3   kHz,    corresponding  to 
an  internal   field  of  0.1    mT.    Therefore  it  is  not  satisfactory  to 
describe  the  high-field  phase  as  a  highly  polarized  paramagnetic  system. 
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As  shown  in  Figure  5.17  the  data  from  the  measurement  at  0.495   T  do 
fit   quite  well  to  a  relation  P  =   PQ   +  AT4  with  A  =  2.66.    The  spin-wave 
velocities  would  be  6.2,   7.8,   or  8.9  cm/sec  for  one,   two  or  three  modes, 
respectively.    In  general,   there  will   be  one  spin-wave  mode  corresponding 
to  each  broken  symmetry. 

5.5     Concluding  Remarks 
The  major  conclusions  of  this  work  are  as  follows: 
1.      For  the  bcc  solid  3He  at  melting  pressure  there  exist  three 
different  magnetic  phases:    paramagnetic  phase,    low-field  phase   (or   u2d2 
phase)   and  high-field  phase.   The  high-field  phase  is  an  ordered  state  as 
indicated  by  the  T     dependence  of  the  melting  pressure. 

It  has   been  suggested  that   the  high-field  phase  has  a  CNAF  (canted 
normal  antiferromagnetic)   structure.   This  work  could  not  give  direct 
information  about  the  magnetic  structure  of  this  phase.    But  the  spin- 
wave  analysis  based  on  the  pressure  measurement  gives   the  spin-wave 
velocities  for   different  modes  which  could  be   useful  for  theoretical 
analysis  of  the  spontaneous  broken  symmetry  of  the  possible  magnetic 
structure  model  for  the  high-field  phase. 

2.      The  magnetic  phase  diagram  of  bcc  solid  3He  at  melting  pressure 
has  been  reestablished. 

Below  0.396   T  the  magnetic  ordering  from  the  paramagnetic  phase  to 
low-field  phase    (u2d2  phase)    is   confirmed  to   be  first  order.    The  entropy 
discontinuity  at  this  transition  decreases   as   the  field  increases. 

From  0.396   T  to  0.410  T  two  phase  transitions  have   been  seen  at  a 
fixed  magnetic  field.    The  transition  between  the  high-field  phase  and 
low-field  phase  happens  at  a  lower  temperature  than  the  transition 
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between  the  paramagnetic  phase  and  high-field  phase.   The  transition 
temperature  between  the  high-field  phase  and  low-field  phase  decreases 
sharply  with  increase  of  magnetic  field.    This  transition  might  be  seen 
at  even  higher  field  up  to  0.4H0  T  according  to  Osheroff's  measurement. 
The  nature  of   this  transition  is  first  order,   as  reported  before.    The 
other  transition  between  the  paramagnetic  phase  and  high-field  phase  is 
not  first  order  from  the  melting  pressure  measurement  of   this  work.    No 
entropy  discontinuity  has   been  seen  at  this  transition  within  the 
sensitivity  of  the  experiment.   The  characteristic  of  dP/dt  shows  a  rapid 
change  around  the  transition  temperature.   The  transition  may  be  socond 
order  or   A   transition. 

The  intersection  of   two  phase  boundaries  has   been  determined  at  B  = 
0.396  T  and  T  =  0.90  mK.    Around  this   point  dB/dT  shows  a  continuous 
change  rather  than  a  discontinuous   change  reported  before. 

The  current  theory  of  multi-exchange  interaction  gives  a 
qualitatively  correct   phase  diagram.    But  there  is  still  a  quantitative 
discrepancy  between  the  experimental  results  and  theoretical  predic- 
tions.   Some  further  experiments  will   be  helpful   to  understand  the  nature 
of  the   different  magnetic  phases   and  their   properties.    Extension  of  the 
experiment  to  measure  the  melting  pressure  of  solid  3He  at  high  magnetic 
field  up  to  10  T  or  even  higher  will   give  more  information  about  the 
high-field  phase  and  verify  the  current  theory  about   the  nature  of  this 
phase.   On  the  other  hand  the  exchange  process  depends  strongly  on  the 
molar   volume  of  solid  3He#    By  studying  the  pr0perties  of  ordered  solid 
JHe  at   different  density  one  will   be  able  to  see  the  density  dependence 
of   the  different  exchange   process.   Mamiya  et  al.    (1981),    Hata  et  al. 
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(1983),  and  Takano  et  al.  (1985)  reported  their  results  related  to  high 
density  solid  3He.  The  experiment  in  hep  phase  of  solid  3He  is  now  going 
on  in  this  lab  to  see  more  closely  the  relationship  between  the  density 
and  exchange  process. 


APPENDIX     A 
THE   RELATION  BETWEEN  SPIN-WAVE  VELOCITY  AND  PRESSURE  MEASUREMENT 

In  analogous  to  lattice  phonons  the  entropy  of  solid  3He  could  be 
written  as 

S3  =  aT3      ,  (A.1) 

where  a  is  a   constant,    T  is  absolute  temperature. 

By  using  the  relation  (A.I),   the  Clausius-Clapeyron  equation 
becomes 


<£,     Ss~Sl     „.«£.  (A2) 

dT  V      -  V,  AV  '  K    'c> 

s  1 


where  Ss,  S^,  Vg  and  V-^  are  the  molar  entropy  and  molar  volume  of  solid 
and  liquid  JHe  respectively,  AV  is  the  molar  volume  difference  between 
solid  and  liquid  JHe.  Here  the  contribution  of  the  entropy  from  liquid 
JHe  below  the  ordering  temperature  is  essentially  zero. 
Integrating  both  sides  of  (A. 2),  we  have 


P  =  P  +  -A-   T4 
0   4AV 


Pn  +  AT4   ,  (A.  3) 
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where  A  =  a/^AV  could  be  obtained  from  the  fit  parameter  of  pressure 
measurement. 

On  the  other  hand  the  specific  heat  per  mole  of  solid  3He  is 


dS 


Compared  with  the  result  of  spin-wave  theory  (Cross  and  Fisher,  1985) 

2   k     1 
2*      kB   T        1 

15  tt°    a  C  J(k)   S 


we  obtain  a  relation  between  spin-wave  velocity  and  parameter  A  from 
pressure  measurement, 


z   <_J >  _  "5  an3 

a     C  J(k)    2ir  kj     s 
a  B 

90  A  AV  ft3  „  .  rs 

=  2-— 14 —  Vs   •  (A. 6) 


where  in  (A. 5)  the  sum  is  for  all  the  spin-wave  modes,  C  is  the  sDin- 

a        K 

wave  velocity  of  the  oth  mode.  The  spin-wave  velocity  will  be  in  cm/sec, 

if  AV  and  Vs  are  in  cm3/mole,  A  in  mbar/(mK)\  and  h  ,  kg  in  cgs  unit. 


APPENDIX     B 
MELTING  PRESSURE  OF  A  FREE-SPIN  SYSTEM 

The  partition  function  of  a  free-spin  system  in  a  magnetic  field  is 
given  by 


m=+J  guRBm 

Z  =     I       exp   (-jpp)      .  (B.I) 

m=-J  B 


where  g  is  the  Lande  g  factor,    pD  the  Bohr  magneton,    B  the  magnetic 

D 

field,   m  the  magnetic  quantum  number,   and  T  the  absolute  temperature. 
Let 

a  =  gUBB/kBT     ,  (B.2) 

then 


m=+J 
Z   =     E       exp    (am) 
m=-J 
=  exp   (-aJ)   +  exp  [-a(J-1)j   +  •••   +  exp   (aJ)      .  (B.3) 


This  is  a  finite  geometric  progression  with  a  ratio  e    ,   the  sum  is 

sinh   ( J+1  /2)a 


Z   = 


sinh  (1 /2)a 
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In  Z  =  in  sinh  (J+1/2)a  -  In  sinh  (1/2)a  .        (B.U) 


Then,  the  entropy  of  the  system  could  be  written  as 


s  .  Nkfl  itttaSL    .  (g>5) 


By  using  the  Clausius-Clapeyron  equation  for  solid  ^He  at  melting 
pressure,   neglecting  the  entropy  of  liquid  ^He,   the  pressure  could  be 
written  as 


P    -   P0    =   /    (f )    dT 


'    (4f-)    dT 


Nk 
fT  ^B      3(T  In  Z) 

J0   AV  3T  a 


Nk 
M"TlnZ 


NkB 

—  T   [in  sinh   (J+1/2)a  -   In  sinh   (1/2)a]      .  (B.6) 


For   3He  J  =   1/2,    so 


Nk 
P  -  PQ   =  Ty-  T   [in  sinh  a  -   In  sinh   (1/2)a]      .  (B.7) 


The  parameter  a  for  3He  coui<j  be  calculated  by  using  the  NMR  data 
v/B   =   32.435   MHz/Tesla    , 
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a  =  iTY-  =  irT =  1-55  t   •  (B-8> 

o  B 


where  B  is  in  Tesla,   and  T  in  mK. 

o 

If  solid  -'He  is  really  a  free-spin  system,   then  the  melting 
pressure  as  a  function  of  temperature  will   be 


NkB 

p  =  po  +  mT  T  tln  sinh  (1-55  t>  "  ln  3inh  (0-775  ?)] 


P  =  PQ  +  AT  ln[exp  (0.775  p    *   exp  (-0.775  f ) ]   ,      (B.9) 


where  A  *  -77— 
AV 


T 
Nk 


APPENDIX     C 
MAGNET  RAMP   PROGRAM 

The   program  is  written  for  HP-41CV  calculator  which  in  turn 
simulates  a  12-bit   programmable  up/down  counter  for  ramping  the  power 
supply  Kepco  AET6-100M  through  Kepco  programmer  SN488-122.    The  rate  of 
magnetization/demagnetization,   step  size,    initial  and  final  step  value 
are  prompted  by  operator. 

There  are  two  operating  modes:    fast  mode(10  V  range)  and  slow  mode 
(1    V  range,   set  flag  01).    The  main  program  RAMP  updates  a  decimal  value 
at  the  selected  rate,   sends  this  value  to  the  output  routine   (HOUT). 
HOUT  converts  the  value  from  RAMP  to  hexadecimal,    formats  the  control 
command,   and  sends   it  to  the  proper   interface  by  a  subroutine  DTOH. 
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APPENDIX  D 
TEMPERATURE  PROGRAM 


10  !  PROGRAfl  -T<Pr>" 

29  FIXED  3 

30  INPUT  "ENTER  PRESSURE  IN  PSI!  ',P 

48  P-P'14.3038          |  CONVERT  PSI  TO  BAR 

38  T«PNN»uton<P,ie,.ei8) 

(8  PRINT  -T  -"ST;"  «K" 

78  GOTO  38 

88  DEF  FNNtyton(Y,xe,D«U«> 

98  Loop:    Y8>FNP(X8> 

»ee  *«<o»<y-y8>/<fnp<x8*ii»  it  •>-¥«> 

lie  x-xe*B«it««R«tio 

128  IF  ABS<R«t1oX.l  THEN  RETURN  X 

138  X8-X 

148  GOTO  Loop 

138  FNENO 

168  DEF  FNP<T> 

178  f  PROG!  PRESS 

168  ,  function:  g^..  th.  m.Uing  pr.„ur.  of  H#_3  b#Joy  th-  m{n ^ 

288  I  REF:  ? 

210  DIN  A<-3:3> 

228  INTEGER  I 

238  I  R»«d  pewvr  s«M«s  coefficient* 

248  READ  *(•) 

Z6B  T2—T#T 

278  T3-T«T2 

288  T4-T2*T2 

298  T5-T«T4 

5?:  p^^flo^>338^*"<■^>r^^^■l:;:::c!^:c^r',<2)#T2+ft<^>•T^*fl<o•T4+"<3>•T5 

328  RETURN  P 

338  FNEND 
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APPENDIX   E 
DATA  ACQUISITION  PROGRAM 


1   !  PROG:  TEMPSfl  >u:  CREC  HRBS  , 

2  |  tunctioh:  p.rror..  ??;.  e.i?.«S5  2Si«..  uh.n  lr^:r;;  ^v^r 

Explorer  III  or  Speei.l  Function  Kr      ■ 
Include*  tint  fro-  clock  interface. 


_   •       ...        ■  -  —  »..•*»  vrifn  (riagrrta  Dv  SRD  from 


a   ,  ?4t*  r»cordtd  by  record  •  (see  H nt  2426). 

9   '  Itll/Vill'l   ^  °i'TPUT  »*•*— "*•  ^«h,r  lh»  PRINT), 

iaa   i  r>.,.         d   (,t  ■*n,P"'«t<on  «ubroutine». 
100   I  Becl»r«t»on.  And  Ini t t ilj 2»t 1 on 

i*  !{jif*i:;::3;:g;rt't,M'«*t»«"«."M«tji«uMci„i,ii„ 

150    M«-"OHX«" 

tee   fixed  a 

f?S   '•J^'"  '*■'•  *nd  »e*'«  functloni- 

f i?    "f.  FN"**n>-INT<MflXce,I^Hnor«).St»p>*K 

f?f    III    ^N"*"<I>-lNT<HlN<L1-,I>'Hnor-).St.p)*IC 

310  C0SU1  Input 

328  FiW-1 

330  COSUB  N*.» 

350  Poi nt •R#cord"0 

360  Printer*© 

500  I  Esttbli.h  .(»e.n«n.ou»  Interrupt. 

310  ON  KEY  t0,3  C0SU1  Start 

520  ON  KEY  tl,5  C0SU1  Input 

330  ON  KEY  #2,7  COSUB  Printer 

540  ON  KEY  .3,?  COSUB  Printer 

550  ON  KEY  .4,7  COSUB  Coiaent 

560  ON  KEY  #5,7  COSUB  Tiae 

570  ON  KEY  #6,7  COSUB  Voltaeter 

380  ON  KEY  07,3  GOTO  Exit 

*?2  'AVtP   E"P1or»p  control  and  tdtt  loop 

610  COSUB  Hold  next 

620  Idle:  URIT  jeo 

630  GOTO  Idle 

650  !  Exit  progr»« 

660  Exit:  OFF  INT  #7 

670  OUTPUT  713;"H2- 

673  STATUS  714;Statu» 

680  DISP  "PROGRAM  BONE" 

690  ENS 

?ii  sr„:hsT«Tusm7M;sr:::r  mm  Exp,or'r  Ho,d'  • "*«• 

"•       OFFN;NT'.J<St4tU,'S>  ™EN  C""-»<— " 
7*9  Point»Point*l 

750      COSUB  Tiae 

?"       Co--•ntS-•■•nt•>>,  ™EN  °UTPUT  Prin*-"!Co-..nt. 
780      GOSUB  Voltaeter 

790  COTO  Expl 

791  Clear_»canner:  STATUS  70»;Statu» 
752  RETURN 
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ENTER  714  USING  Nora_f  era;  Vbi  t  <  I  >,Hbi  t  CI>,Hstep<  I ) ,  V0<  I  ),H8C  I  >,  Vnora 


899   I  Enter  endpoints  of  intervals 

■!"e!>lnUt:INPUT  "ENT"  ,ECIMN,NC  "ND  EHB  OF  BASELINE  AVERAGING  INTERVAL  Cu,.t>: 

•  3«        nntlln-ll"  ,ECIHNINC  "ND  ENB  0f    INTEGRATION  INTERVAL  Cu..c  >:  ',  II,  In 

848       RETURN 

"5   s!Jr""uTrCTr7?§rH2?    <nt'»r*tlon   *   '•«•"■   t^or.r   t.   Held 

ere      return 

968   I  Printer  definitions  and  coaaent  Input 

910  Printer:  Printtr-lSeNOT  Printer 

"*0  IF  Printer-8  THEN  DISP  "OUTPUT  TO  THERMAL  PRINTER" 

939  IF  Printer-16  THEN  DISP  "OUTPUT  TO  CRT" 

946  RETURN 

956  Coaaent:  EDIT  "WHAT  IS  COMMENT  ?", Coaaent* 

968  RETURN 

1868  !  Determine  how  auch  aeaory  is  accessed 

ieie  Expi:  output  7i5;*D4" 

1913       ENTER  714  USING  *N" ; Track .Memory 

1828        RESET  714 

1838        0trsaINT(H«Bory4>*l 

10*0        Mt»loc»Me*ory  NOD  4 

1838        IF  Heaory<6  THEN  Six* 

1868  Qtrs-4 

1078  Healoc'6 

1080  Size:  Q"2»Qtrs-l 

1898       Scale«4'0trs 

1180  !  Read  nonillittion  Information 

1110  Norm: DIM  Vbi t <7> , Hbi t <7> , Mst tp(7>, V8C7> , H8<7> , Vnora<7> ,Hnora<7> , St ep<7> 

1128  Nora_fora:  IMAGE  3(N> , 2C3N> , 2<7N7 

1138  OUTPUT    713j"N2* 

1148  FOR    1-8   TO   0 

use 

CD.Hnon 

1168  NEXT    I 

1178      RESET  714 

nee  hat  step-M»tep'<scaie> 

1198  MAT  Hnora"Hnora»(lE6> 

1288  !  Read  digitized  signal 

1218  Read_data:REDIM  DataC 1624eQtrs-l> 

1228  OUTPUT  713; *D2" 

1239  ENTER  714  BFHS  2048*Qtr»  USINC  "t, Y" ; Dat a<«> 

1248  RESET  714 

1308  !  Separate  data  as  recorded  and  Integrate 

1318  Sep:  FOR  K«8  TO  0 

132B  IF  Step(K)-e  THEN  Hextk 

1330  Step«Step<K> 

,349  Seg*-Memloc*49*KeCl*(<Qtrs-2>  AND  (Step>2>>> 

1359  FOR  J-K  TO  0  STEP  HAX<Step,i> 
1353  IF  StepCJ>OStep  THEN  Data  »rr 

1360  Step(J)«8  ~ 
«3?e          NEXT  J 

1388  COSUI  Liaita 

1398  Nextk:  NEXT  K 

14ee  !  Enable  interrupt,  put  Explorer  in  Live  aode,  and  return. 

1418  Ho1d_next:  ON  INT  87,3  GOSUB  Srq 

1*20  CONTROL  MASK  7^128 

1438  OUTPUT  713; 'HB' 

l«<0  RESET  713 

1*50  CARD  ENABLE  7 

1460  RETURN 

1388  I  Signals  recorded  in  different  size  aeaory  groups  (e.g..  H2  1  04) 

1318  Data_err: I 

%V1>  ~0UTPUT  Winter  USING  -VK";"DATA  INCONSISTENT  FOR  •lHSC4'Hst ep;  I  UCHR»< 

1338      GOTO  Nextk 

2888  I  Set  up  interval  endpoints 

till   U'U»SJ-1-4""""«»<'<>-»  I  •  of  point.  1n  aeaory  group 

2628        HnoraaHnora(K) 
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««•       if  n«<.:"«N  n:::Lr:.ir.9:;;4on  •■•• *-■  "  *«•• « « 

«=•       n.-rNHLcii; to  ,,n9th  fp°"  egrr*nt  •»••*•  * r  .<;*!;. 

ff!2  Ih1-FNH(n<I§*nnt>  !  Th***    Hn"   •"»«""•   »h«t 

«■•■  «to-FNP1«K<«i>  •ndpoint*    art    vtlhtn 

"'2  «M.FNH(n<-,n>  '        r*n»»   of   «*•*•• 

,!:!  IF   "•<••   THEN   ■lo-Ilo  ,  t,    ,      , 

llll    .    «  IF   Rn<"M   ™EN   "h«-   hi  ^l^   !    •W*P*8,n«   •"dpoint. 

2   M    !    »„„.,„,   „tp   „„   „,„  t  •         T.    int.gr.tion  .„dpo(nt,. 

llll  "   »*•"•>•>    THEN   H..d.r 

2»s»  St, p.) 

2J7I  m..-      c  PM"  <PftRT"»-  »«Tfl>  :■ 

till    t4d,r-^?^t*"",t0^»^.P;ntCHR,(S,g.> 

ss  ^•;£:"- -™~  2:1:="*  8""E ••"• 

2248  NEXT    I 

2m  enn'SW  •b,°'ut'  U4,u*  of  •'•"•' 

"§S  ro?J*^!  T0  n°  STEP  -*"p  '" 

l"l  Ch,.ch„MScD.t.<n-RW.,P 

Ilirv  f™^Bt."r  INTECRflL  or  sicNflL  FRon  ■«'-"„.>,..,..«  ro  .;FNTah1 

2«J!  ^x\.vnd:r  —  -*«-  *- «...  ...0^  oroup 

2«2   RETURN1  "  •■•"•"-!^««.  S t«,Ti..,|ChU,Volt.<.> 

*■••  1  R,*d  «•*•  *nd  M.. 

3818  Ti»t:  OUTPUT  6; "R» 

3928        ENTER  SjNonth, D*y, T..,» 

■TinE..T<..rPUT  *"""  USINC  T*-^-a.;-I»«TE.,f,o„th,...,D4v)>.polNT  ,.fpejnt 

3»ee  !  Scan  !oek-fn» 

31  16  Volmtirj  OUTPUT  789: 'DCVS-*. 

3}2|  ENTER  r.";.?"?.;  '  T«"  3421ft  to  r..d  «h.  2.5 

.XT'^'T^l^  —  —  -*.. 

is  Z.  ^P-":" »"«?:-"'"  ""■•'"•• 

«i«2    ■    -,  '!   ?ok"a    ™«   RETURN 


4, an     1    m  .    "-«-■     men    RETURN 

4,"      Prop#r.!'!#  «••■  "O.  •*<« 


,,3B  R*ply>»«>y» 

«>«•  INPUT  -SHOULB  I  CREATE  THE  FILE  V.R.ply, 


126 


4,s»  IF   UPCt<Rtp1y»tl;tlX>"Y"    THEN   H*a* 

4»<»  CREATE   FiU6f:T14a,23S,48 

4»7«  COTO  A»»ian 

4266  I  List  rtl*  tiblt 

4218  F1U_t«bJ»:  PRINT  LIN<1>,»FILE  6", 'NRHE'.LINC  1  > 

«228  FOR  I-J  TO  16 

*>23t  PRINT  l,Filit«<I) 

*>248  NEXT  I 

•236  RETURN 

5686  I  Ffti  atnipul ttlon  sain  routine 

3618  Ntnipu) *tt:   I 

S188  I  Loe*tt  Input  flit 

sue  Fiittn:  fi  i*-2 

3126  COSUI  N*at 

3136  Vn-4 

5148  INPUT  -HOW  HANY  VOLTAGE  READINGS  ARE  IN  EACH  RECORD  ?".Vn 

5136  REDIH  V<Vn-l> 

5166  R«ply«»"Y" 

5179  INPUT  "IS  THIS  A  RAM  DATA  PILE  ?  <  Y  •  N  >",Rtply« 

3>86  IF  UPC«(R»ply»Cl;i:>-"Y"  THEN  R*u 

5206  !  R»«d  procttttd  file 

3218  ON  END  62  COTO  Flit  tnd 

5226  FOR  W  TO  1686 

3236  READ  62,Rtc jPt < I >,T», Rtlt < I ),Chi < I ), V<») 

3248  NEXT  I 

5258  COTO  Ft  1*  tnd 

5388  '    R»»d  r*u   dut'riU 

3318  R*y.    READ    82, 1 ;Pt < 1 > , S», T*, Chi < 1 > , V<«> 

5328  L»*tti»t-Ti«te-(VAL<T*tn)«66*VAL<T*t4J>>#6e*VAL<T»t7]> 

5336  ReUCl>"8 

5338  ON   END   82   COTO   Fit*   tnd 

3368  FOR    1-1    TO   236 

3378  READ  62,  Rtc  ;Pt  <  I  >,  S8,T8,  CM  <  I  >,  VCO 

3388  T1»«»(VALCT«Cn>»6e*VAL(T«C43)>«68*VALCT«C73) 

3398  IF  Ltstt  inOTtat  THEN  Ti  ««8«Ti  ■•6-86488 

3488  RtHU>«(Ti»f-Ti«»e>/36e0 

3418  l»»t i i»««Ti«» 

3428  NEXT  I 

5438  Filt_tnd:  L»»t"l-1 

3468  RETURN 


APPENDIX     F 
ENTROPY  PROGRAM 


10  !    PROCRflH    "SPLOT"    reads    *   data  ftlt,    calculate*    tht    entropy   of    solid  He3 

,-  :                                             *nd   P'ots    the    prrtiurt    vi.    tin*,    pritturt   vs.    Umpirttun, 

ZZ  *nd   tntropy    of    solid    He3    vs.    temperature. 

48  ■    Variable   table: 

50  I    Ti-tiat 

68  !  T-te«iperature  from  the  fit  of  NMR  vs.  tine 

78  !  Tc-temperature  of  bottom  strain  gauge 

80  !  Chi-NHR  temperature, 

98  I  Er-difference  between  Te  and  Chi 

188  !  S-entropy  of  solid  He3, 

118  !  H-speeific  heat  of  solid  He3 

120  !  C-pressure  of  Me3  in  the  sample  cell  in  Bar 

130    OPTION  1P.SE  1 

140  L-450 

150  Sin  PlottertCCJ 

lit  I1m*l    T'(*Se>>T<45e>»Te<45e>.ChiC450),Ert450),S(450)lSl<450),H<450) 

170  DIN  CC430) 

180  INPUT  "Enter  name  for  data  file", Filet 

190  ASSIGN  01  TO  Pi  let 

288  ON  END  tl  GOTO  250 

210  FOP  I«J  TO  L 

220  READ  «1, I;Ti<I>,C(I),Chi<I),Te<I>,ErCn 

238  !  PRINT  USING  "3D, 2x, 3D. 2D, 2X, 30. 3D" : I , Ti < I > . CC I > 

240  NEXT  I 

230  n-I-1 

260  PRINT  "r>",M 

270  FOR  !■!  TO  H 

280  C<l>-<C<I)-34. 3380). 1000    .  P-P.  i„  „..,.,   p4  ,,  th,  prt„ur,  ,t  « 

298  NEXT  I                       trasition  based  on  Creywal 1 's  data 

308  !  Change  t i m*  to  temperature  according  to  the  fit  of  NfIR  temperature  us. 

t  i  me. 

318  INPUT  "Enter  the  coefficients  of  the  fit  f unet i on", R0.P1 

320  FOR  I"l  TO  H 

338  TCI)"P.e*P.l«Ti<I> 

340  NEXT  I 

350  !  Routine  for  smoothing  the  pressure 

368  Imax-n-l 

378  C1«C<1) 

388  FOR  1-1  TO  Imax 

390  Cnew"<Cl*2*C<I>*C<I*l>>».23 

400  Cl-CCI) 

410  C(I)-Cnew 

420  NEXT  I 

430  C<n)-CC10*C<n>>*.2S 

440  FOR  I"l  TO  H-« 

450  S<I)--<C<I*«)-CCI>)x<T(I*6)-T(I)).. 13^3.763   I  Constant  .13   -3.763 

puts  S  in  unit  Rln2  for  one  mole  He3 

460  NEXT  I 

470  M-H-6 

460  !  Add  a  small  contribution  of  liquid  He3  to  the  entropy 

490  G-6.636E-3     !  C'Camma* 10E-3- 1 n2,   Camma-4.6 

See  INPUT  "Give  the  point  number  of  P,  t  ransi  t  i  on"  ,  N 

III  «??':?,I°  "   '  "dd  S'  T0  S  t0  ■•«  s»  •*•«•  *"•  "  transition 

5*0  S<n«S<I>*C4(T(I>*T<I*6)>'2 

530  NEXT  I 
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548    FOR  1-1  TO  H 

550    ll(I).(S<l^).$(I»/(T(I»«M<I))    I  Ctleultti  dS»'dT.R.ln2 
III        H<I>'SH!).(T(I)»T(I»«))/2    |  Spteifie  htat  <H>  for  solid  Hc3 

see   next"!  USINC  "3B,3x,3D-41I'3><''131,-3D'3x.2D-31,":i.t<»>.c<:).s<i> 

590    Taax»T<l> 
608    TitniT(l) 

6ie   c»*x«co> 

620    Ciin>C(l) 

638    5n».$(|) 

6«e   s»in«scj> 

650    FOR  1-1  TO  H 

66e    IF  TCI>>Taax  THEN  Taax-Td) 

676         IF   TdXT.in   THEN   Tain»T<I> 

688         IF   CCI>>Caax    THEN   Caax-CU) 

698         IF    CUXCain    THEN   Cain«C<I> 

788    IF  SCI>>Saax  THEN  Saax-S(I) 

716         IF    SdKSiin    THEN   Sain-S(I) 

728    NEXT  I 

738    PRINT  •TBin"';T«in,"T«*x»-;T««x 

748    PRINT  "C«in-"JC»in,"C»»x-"jC«i*x 

758    PRINT  "Sain»";S«iin, "Smax-" ; Saax 

768    INPUT  "GIVE  THE  NAME  YOU  USE  FOR  PLOTTING  (CRT'9672A>  •  ,  PI  ot  t  »r« 

778    IF  Plott»r$«"CRT"  THEN  798 

768    IF  P1ott«r»-"9872A"  THEN  828 

798  PLOTTER  IS  13, "GRAPHICS" 

888  GRAPHICS 

616  GOTO  678 

628  PLOTTER  IS  "9872R" 

638  OUTPUT  785; "VS5" 

648  LIMIT  8,279.4,8,215.9  "For9872fl 

858  LOCATE  14.5,128,12,94  iFor  9872R- 

668  GOTO  698 

678  LIMIT  8,184.47,8,149.82  «For  CRT 

688  LOCATE  14.5,128,12,98  !For  CRT 

III  aatul  lZiV'    th*  w*,u»*  of  X"*".X»ax, Vain, Yaax" , Xai n, Xaax, Vain, Yaax 

908  PRINT  Xain.Xaax, Yain, Yaax 

918  INPUT  "Ciut  th.  valutt  Of  X»t»p,Y»t»p",X»t»p,Y»t.p 

928  PRINT  X»ttp,Y»t«p 

938  SCALE  Xain.Xaax, Yain, Yaax 

948  AXES  Xst»p,Y»ttp,Xain,Yain 

958  AXES  X*tap,Y»ttp, Xaax, Yaax 

9c8  FRAME 

.  978  COSUB  Lxaxtft 

968  GOSUB  Lyaxts 

998  LORG  5 

1888  MOVE  <Xwax*Xain>'2,Yain-<Yaax-Y»in>'ll 
1818    LABEL  USING  "K" ; "T.aperat ur»<aK> " 

1828  !  LABEL  USING  "K" ;  "Ti  M(hn)  • 

1838  MOVE  Xaax-(Xaax-Xain>^4,Yaax-<Yaax-YBln)»3^4 

1848  INPUT  "  Uhat  it  the  aagnvtic  fitld'".F 

1858  LABEL  USING  "K";"B  «";F;"T»»1a" 

1868  MOVE  Xain-(Xaax-XBin>^e,<Yaax+YBin)^2 

1878  DEC 

1888  LORG  5 

1898  LDIR  90 

1188  I  LABEL  USING  "K*;"P-Pa  (aBar)" 
1110    LABEL  USINC  "K";-S»'R*1n2" 

1120  I  LABEL  USING  "K"; "dS»'dT»R»ln2" 

1130  I  LABEL  USINC  "K" ; "Cv'R«ln2" 

1140  FOR  la]  TO  N 

1150  I  PLOT  T1(I>,C<I) 

1160  I  PLOT  T<I>,CU> 

1170  PLOT  <T(I>*T<I+6)>'2,S<I> 

1160  I  PLOT  <T<I>*T<I*C>>'2,SKI> 

1196  I  PLOT  <T<I)*T<I*6>>^2,H<I> 
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1269  PENUP 
121*   NEXT  I 
1220   END 

1230  Lx*xt»:   I 
1240  CSIZE  3 
12S0  LORC  S 
1260  LDIR  0 

1270  FOR  Xpo»iHon-X«in  TO  Xa*x  STEP  X»ttp 
1260  HOVE  Xpo»U1on,Y«tn-<y««K-Y>1n>x22 
1290  LflJEL  USING  "B. »■; XpoH t ion 

1300  NEXT  Xpositton 
1310  RETURN 

1320  Ly*xc»:    I 
1330  CSIZE  3 
1340  LOIR  0 
1350  LORC  0 

1369  FOR  Ypo»it ion-Yain  TO  Y«»x  STEP  Y»tep 

1370  HOVE  Xnin, Ypotitien 

1380  LABEL  USINC  T13D.DX";Ypo»ition 
1399  NEXT  Ypotitten 
1499  RETURN 
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